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Abstract 


The  overall  goal  of  this  work  is  to  develop,  implement,  and  evaluate  novel  three-dimensional 
acquisition  orbits  used  in  conjunction  with  a  dedicated  emission  mammotomography  system,  to 
be  used  as  a  potential  secondary  diagnostic  tool.  This  novel  nuclear  medicine  based 
molecular  imaging  technique  utilizes  camera  trajectories  that  sample  a  pendant,  uncompressed 
■breast,  which  are  evaluated  with  respect  to  various  breast  shapes  and  sizes.  Phantom 
]  measurements  have  demonstrated  that  lesions  near  the  chest  wall  can  be  visualized  for  both 
large  and  small  sized  breasts,  and  furthermore,  that  axillary  imaging  is  possible  even  with 
minimal  system  tilt.  The  use  of  additional  bed  shielding,  simulated  here  with  a  lead  sheet, 
was  also  investigated.  Preliminary  results  suggest  the  lead  shield  may  hinder  lesion 
detectability.  Future  work  includes  the  addition  of  a  patient  bed,  which  is  expected  to 
further  complicate  positioning,  and  which  may  or  may  not  include  lead  lining  pending  further 
;  studies .  _ _ _ _ 
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Introduction 


The  overall  goal  of  this  work  is  to  develop,  implement,  and  evaluate  novel  three-dimensional  acquisition  orbits  used  in 
conjunction  with  a  dedicated  emission  mammotomography  system,  which  may  allow  for  earlier  detection  of  a  primary 
breast  cancer.  The  tomographic  gantry  central  to  the  ASET  allows  for  camera  positioning  anywhere  about  an 
uncompressed,  pendant  breast  in  inverted  hemispherical  volumes  [1,  2],  The  key  benefits  of  imaging  in  this  geometry 
include  the  increased  separation  of  the  breast  tissue  from  the  nearby  body  containing  background  and  the  ability  to 
image  an  uncompressed  breast.  One  main  advantage  of  3D  orbits  that  utilize  increased  polar  tilt  is  then  an  increased 
viewable  breast  volume,  potentially  into  the  chest  wall  region  [1],  However,  with  increased  tilt  angle,  the  camera  may 
directly  view  cardiac  and  hepatic  regions,  which  could  contribute  primary  background  contamination  to  the  2D 
projection  images  of  the  breast.  Scintimammography  phantom  imaging  studies  indicate  that  scatter  from  cardiac  and 
hepatic  background  can  significantly  contaminate  the  primary  lesion  signal  [3].  Studies  with  clinical  SPECT  gamma 
cameras  using  simple  tilted-parallel-beam  (TPB)  vertical-axis-of-rotation  (VAOR)  orbits  [4],  and  simulations  of 
somewhat  more  sophisticated  orbits  [5]  have  shown  the  potential  effectiveness  of  dedicated  SPECT  breast  imaging 
approaches  that  use  non-standard  camera  orientations.  Lesion  contrast  and  SNR  improve  by  nearly  a  factor  of  three 
when  using  a  clinical  camera  with  tilted-parallel-beam  as  compared  to  low  noise  planar  (2D)  images  [4],  However, 
current  clinical  SPECT  scanners  are  not  only  limited  by  the  ROR  that  they  can  achieve,  but  also,  by  the  trajectories 
possible  with  the  gantry.  Through  its  combined  fully  computed  controlled  variable  radius-of-rotation,  polar  and 
azimuthal  angular  positioning,  the  ASET  has  opened  up  a  new  class  of  volumetric  orbits  for  tomographic  acquisition 
not  previously  possible  with  clinical  scanners. 

Body 

Task  2.  Development  and  investigation  of  novel  orbits.  (Please  see  Appendix  A  for  original  statement  of  work). 


Fig.  1.  (LEFT)  Within  each  box,  the  left  plot  shows  polar  camera  tilt  as  a  function  of  the  azimuthal  angle  for  each 
examined  orbit;  the  right  plot  shows  the  radius  of  rotation.  (RIGHT)  Voxelized  Orlov  Volumes  of  the  orbits.  VAOR 
and  PROJSINE  0-45°  are  shown  at  the  95%  threshold.  Because  TPB  45°  and  PROJSINE  15-45°  do  not  meet  Orlov's 
criteria  and  thus  do  not  completely  sample  the  breast,  they  are  instead  shown  at  their  maximum  threshold.  The 
threshold  is  the  percent  solid  angle  that  is  sampled. 

Volumetric  acquisition  orbits  that  more-nearly-completely  sample  the  pendant,  uncompressed  breast  were 
investigated,  with  respect  to  various  breast  shapes  and  sizes  (as  a  continuation  of  Task  2(a)).  Several  orbits  were 
implemented  to  compensate  for  the  positioning  demands  that  are  required  for  different  breast  sizes  (Fig.  1);  while  a 
general  orbit  design  may  remain  the  same  between  two  different  breasts,  the  absolute  polar  tilt  range  and  ROR  range 
may  vary.  Thus,  these  orbits  were  designed  with  a  specific  focus  on  orbit  parameters  of:  viewable  breast  volume, 
radius-of-rotation,  location  of  the  center-of-rotation,  and  polar  tilt  angle,  considering  calculated  Orlov  Volumes  [6]. 
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Calculated  Orlov  Volumes  are  shown  in  Fig.  1.  We  evaluated  this  newer  set  of  developed  orbits  with  measurements  of 
lesion  contrasts  and  signal-to-noise  ratios  in  anthropomorphic  breast  phantoms  with  and  without  additional  torso 
backgrounds,  and  (2)  examining  the  effect  of  various  concentration  ratios  and  various  lesion  sizes  on  image  quality,  to 
characterize  the  limitations  of  the  system's  performance.  For  detailed  methods  and  results  from  the  breast-only 
studies,  and  preliminary  reconstructions  of  the  bed  shielding  studies,  please  see  Appendix  B. 

Detailed  descriptions  of  the  methods  for  the  additional  bed  shielding  experiments  are  given  in  Appendix  B. 
Photographs  of  the  setup  are  shown  in  Fig.  2.  Newly  reconstructed  results  of  the  bed  shielding  studies  are  shown  in  Fig. 
3. 


Fig.  2.  Photographs  of:  (LEFT)  1060  ml  breast  phantom  attached  to  the  torso,  with  lesion  sizes  as  labeled;  (RIGHT) 
Two  vantages  of  the  lead  sheet  used  to  shield  the  'patient'. 


PRO J SINE,  PROJSINE.  Without  lead  sheet  With  lead  sheet 


Fig.  3.  Images  from  reconstructed  1060  ml  breast  +  filled  torso  data  (OSEM,  2.5  mm3  voxels,  2nd  iteration,  calculated 
attenuation  correction,  Hann  filter  applied)  with  and  without  the  lead  sheet.  (LEFT)  3-pixel  wide  profiles  drawn 
through  the  2.27  ml  lesion.  The  lesion  near  the  chest  wall  is  not  visible  with  VAOR  (dashed  white  arrow)  and 
artifacts  move  further  into  the  breast  for  the  PROJSINE  15-45°  case  as  compared  to  PROJSINE  0-45°  (smaller  white 
arrows).  (RIGHT)  More  background  activity  is  visible  in  the  sagittal  views  without  the  lead  sheet;  increased  tilt 
views  the  heart  and  liver  more  directly. 
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Employing  3D  data  acquisition  camera  trajectories  with  increased  polar  camera  tilt,  we  have  demonstrated  that 
lesions  near  the  chest  wall  can  be  visualized  for  both  large  and  small  sized  breasts  (Fig.  3  and  Appendix  B).  As 
expected,  the  shapes  of  the  reduced-artifact  regions  in  the  reconstructed  images  correspond  to  the  respective  shapes  of 
the  calculated  Orlov  Volumes  (Fig.  1).  Because  previous  studies  have  suggested  the  tomographic  sampling  criteria  for 
completeness  may  not  be  as  strict  as  initially  anticipated,  we  relaxed  the  restrictions  on  the  camera  polar  tilt.  This 
allowed  the  COR  to  move  further  into  the  breast,  which  increased  the  reduced-artifact  region  of  the  breast  (Fig.  3), 
thus  allowing  easier  visualization  of  lesions  near  the  chest  wall.  If  a  lesion  is  suspected  near  the  chest  wall,  then  to 
confirm  diagnosis  of  a  lesion  in  that  region  one  may  utilize  more  camera  tilt  and  relax  the  sampling  criteria  (e.g., 
PROJSINE  15-45°);  to  determine  whether  a  smaller  lesion  and/or  lesion  of  low  radiopharmaceutical  uptake  is  present 
anywhere  in  the  breast,  an  orbit  such  as  TPB,  with  its  decreased  ROR,  may  be  more  appropriate. 


• — 

--  VAOR.  without  lead  sheet 

— 

~  TPB  45,  without  lead  sheet 

— 

-  Projsine  0-45.  without  lead  sheet 

— -A— 

-  Projsine  15-45,  without  lead  sheet 

. O . 

VAOR.  with  lead  sheet 

. 0 

TPB  45.  with  lead  sheet 

— © . 

Projsine  0-45,  with  lead  sheet 

. -A- 

Projsine  15-45.  with  lead  sheet 

Fig.  4  Contrast  versus  SNR  for  the  reconstructed  breast  +  torso  data  and  the  2.27  ml  lesion.  Iterations  labeled  next  to 
VAOR.  Patient  shielding  yields  slightly  lower  SNR  and  contrast  values. 

Preliminary  results  suggest  that  the  lead  shield  may  actually  hinder  lesion  detectability,  as  the  SNRs  and  contrasts 
slightly  decrease  (Fig.  4).  One  possible  cause  of  this  decrease  may  be  from  scattered  counts  of  photons  that  emerge  from 
the  breast  and  get  deflected  back  to  the  detector  by  the  shield.  This  trend  may  or  may  not  be  indicative  of  the  true 
nature  of  the  effect  of  the  shield;  with  only  one  measurement  of  each  case,  the  current  measurements  are  likely  within 
error  bars  of  each  other.  Furthermore,  these  measurements  do  not  account  for  the  contralateral  breast,  which  would  also 
contribute  to  the  background.  The  lead  shield  may  not  be  necessary,  given  a  judicious  choice  of  the  acquisition  orbit  and 
subsequent  design  parameters  to  minimize  direct  views  of  the  heart  and  liver,  and  utilizing  the  superior  energy 
resolution  of  the  CZT  camera  to  diminish  the  effect  of  secondary  scattered  background  contamination.  Because  a  lead 
shield  may  preclude  close  proximity  positioning  of  the  detector,  however,  this  preliminary  result  will  be  further 
studied,  with  the  phantoms  positioned  on  a  patient  bed  rather  than  freely  suspended  from  a  phantom  holder. 

In  order  to  determine  which  acquisition  and  reconstruction  parameters  may  be  most  optimal  (Task  2(c)),  software  was 
developed  to  process  previously  acquired  list-mode  data.  One  goal  of  the  post-processing  was  to  examine  the  resulting 
normalized  photopeak  counts  extracted  from  the  summed  energy  spectra  as  a  function  of  the  azimuthal  position  (Fig.  5). 
The  increased  photopeak  counts  from  45<0<135  (where  0  is  the  azimuthal  angle  of  the  rotation  stage)  are  due  to 
contamination  from  the  heart  and  liver,  which  provide  both  primary  and  secondary  scattered  counts.  Increased  polar 
tilt  angles,  which  view  the  heart  and  liver  more  directly,  are  a  direct  cause  of  the  increase  in  the  number  of  photopeak 
counts.  The  larger  polar  tilt  values  correspond  to  the  increases  in  photopeak  counts  (compare  Figs.  1  and  5,  TPB,  45°, 
PROJSINE  0-45°  and  PROJSINE  15-45°) 
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Fig.  5.  Normalized  photopeak  counts  as  a  function  of  the  azimuthal  position,  which  includes  all  valid  events  from  the 
original  list-mode  data  files  within  a  ±8%  energy  window. 

Additionally,  we  ran  studies  with  simulated  repeated  measurements  to  determine  the  variability  in  SNR  and  contrast 
for  each  orbit  given  in  Fig.  1,  given  various  reconstruction  and  acquisition  parameters.  Initial  reconstructions  typically 
do  not  include  calculated  attenuation  compensation  nor  additional  post-processing  with  3D  filtering;  subsequent  re¬ 
reconstructions  included  these  components,  adding  each  independently  (i.e.,  first  add  attenuation  compensation  in  the 
routine,  then  apply  3D  filtering).  Results  from  the  different  reconstruction  techniques  are  given  in  Fig.  6  and  suggest 
that  for  orbits  that  employ  increased  polar  tilt,  contrasts  are  better  without  attenuation  compensation,  while  SNRs 
remain  about  the  same.  In  all  cases,  as  expected,  filtering  increases  the  SNRs  and  decreases  the  contrasts. 


FIG.  6.  (LEFT)  Images  from  reconstructed  results  using  various  reconstruction  techniques  (all  shown  at  the  2nd  iteration). 
For  the  Plain  and  Attenuation  only  images,  there  are  3  summed  slices  to  reduce  noise.  The  filtered  images  do  not  have 
summed  slices.  Profiles  are  9  pixels  wide.  The  dotted  line  represents  the  profile  center  (4  pixels  on  either  side).  Plain  = 
no  attenuation  correction  nor  Harm  Filtering.  AC=Attenuation  Correction.  Harm  filter  applied  at  0.7*Nyquist. 
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FIG.  7.  (LEFT)  Contrasts  as  a  function  of  SNR  for  the  PROJSINE  0-45°  orbit,  using  the  full  counts  as  acquired  during  the 
scan,  and  then  re-processed  with  only  half  of  the  counts. 


Subsequent  studies  included  varying  the  number  of  projections  and  number  of  counts  of  the  PROJSINE  0-45°  orbit,  to 
determine  whether  the  angular  sampling  plays  a  key  role.  Results  from  these  experiments  are  shown  in  Fig.  7.  "Full 
counts"  utilizes  the  full  number  of  counts  acquired  during  the  scan;  "half  counts"  includes  projection  files  post-processed 
to  reduce  the  total  number  of  counts  by  half.  To  create  a  projection  file  with  half  the  number  of  projections,  every  other 
projection  angle  was  written  out  to  a  new  file.  The  "half  counts"  reconstructions  provide  lower  SNRs  and  contrasts  than 
their  "full  counts"  counterparts,  given  equal  numbers  of  projections.  The  comparison  we  are  most  interested  in  is  of  the 
150  projections,  half  counts  to  75  projections,  full  counts  (equivalent  total  number  of  counts).  The  150  projections,  half 
counts  result  in  lower  SNRs  and  contrasts  than  the  75  projections,  full  counts  setup.  However,  at  increasing  iterations, 
the  SNRs  and  contrasts  are  within  error  bars  of  each  other.  Because  of  the  statistical  reconstruction  algorithm,  the 
angular  sampling  may  not  play  as  large  a  role  as  initially  anticipated. 


Fig.  8.  (LEFT)  Two  vantages  of  the  positioning  problem  that  arises  from  the  patient's  arm.  The  camera  must  be  at  an 
appropriate  polar  tilt  angle  and  ROR  to  get  around  the  arm,  which  corresponds  to  an  azimuthal  location  of  0=~3O°. 
Arrows  point  to  limit  of  polar  tilt  and  ROR  at  this  azimuthal  location.  (RIGHT)  Photographs  of  various  degrees  of 
system  tilt,  (left)  0°  system  tilt  (relative  to  the  horizon;  solid  lines).  The  torso  is  tilted  -30°  (dashed  line)  so  that  the 
nipple-chest  axis  is  perpendicular  to  the  horizon,  (right)  15°  system  tilt.  The  torso  itself  was  tilted  an  additional  15° 
(dashed  line),  moving  the  COR  into  the  outer  quadrants  of  the  breast. 

In  another  set  of  measurements,  two  system  tilt  angles  were  examined  to  determine  whether  axillary  imaging  was 
achievable  in  this  configuration.  With  the  additional  complications  that  arise  with  patient  positioning,  we  also 
wanted  to  determine  whether  system  tilt  would  significantly  alter  acquisition  orbit  design.  Data  was  initially 
acquired  about  suspended  anthropomorphic  breast  and  torso  phantoms,  without  a  true  patient  support  palette.  An 
example  of  the  positioning  problem  that  arises  from  the  patient's  arm  is  shown  in  Fig.  8.  With  multiple  degrees  of 
freedom,  the  parameters  we  may  manipulate  to  get  around  the  arm  include  camera  tilt  and  the  ROR.  For  this 
particular  polar  tilt,  the  camera  was  as  close  as  it  could  get  without  hitting  the  'patient'.  The  options  are  thus  (1)  to 
maintain  the  polar  tilt  angle  but  significantly  increase  the  ROR  to  move  the  camera  completely  away  from  the  arm; 
or,  (2)  maintain  a  small  ROR  by  instead  increasing  the  polar  tilt.  We  chose  to  increasing  the  polar  tilt,  resulting  in 
more  direct  vantages  of  the  heart  and  liver;  however,  spatial  resolution  degradation  was  minimized.  In  order  to 
maximize  the  viewable  breast  volume,  the  breast  was  lowered  into  the  FOV  such  that  the  COR  was  as  far  into  the 
breast  as  practicable.  The  only  rotational  motion  was  of  the  torso  itself,  with  a  15°  rotation.  Thus,  in  the  case  of 
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increased  system  tilt,  the  COR  was  actually  located  in  the  left  half  of  the  breast  (Fig.  8,  right),  which  we  predicted 
might  allow  for  axillary  imaging.  Three  130  |iL  balloons  (Harvard  Apparatus,  Holliston,  MA)  (Fig.  9,  upper  left)  were 
used  to  simulate  uptake  in  the  axillae.  Each  was  filled  with  ~14  pCi  of  Tc-99m  and  subsequently  placed  in  wells  ~2.6  cm 
deep  (Fig.  9,  right).  In  addition,  a  2.3  mL  lesion  (Radiology  Support  Devices,  Newport  Beach,  CA)  was  placed  near  the 
center  of  the  935  mL  breast  phantom,  anterior  to  a  large  "cold"  plug.  Projection  data  were  acquired  for  20  minutes  using 
both  simple  and  complex  3D  orbits;  specific  parameters  for  each  of  four  orbits  are  shown  in  Fig.  9. 


System 
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of 
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Fig.  9.  (LEFT)  Photographs  of  the  phantoms,  including:  (upper  left)  Three  130  pL  balloons  to  simulate  uptake  in  the 
axillae  and  one  2.3  mL  lesion;  (lower  left)  935  mL  anthropomorphic  breast  phantom  With  the  2.3  mL  lesion  located  near 
the  breast  center,  anterior  to  a  large  "cold"  plug;  (right)  Location  of  the  three  simulated  lymph  nodes  (breast  has  been 
removed).  Wells  are  ~2.6  cm  deep  at  each  numbered  location.  (RIGHT)  Parameters  Used  for  each  Breast-Only 
Acquisition,  with  all  orbits  acquired  over  a  360°  azimuthal  range. 


Sys.m.TUt.=ffl>  S.v.3fl!i.Tili  =  l?‘ 


PRO.lSINh, 

TPB,  45*  15-45*'  TPB,  60“  COMPLEX 


Fig.  10.  Images  from  OSEM  reconstructed  data  (8  subsets,  3  summed  slices,  2nd  iteration).  2-pixel  wide  profile  drawn 
through  the  2.3  ml  lesion.  Bottom  rows  are  MIP  images  through  the  breast  volume,  to  aid  in  visualization  of  the 
axillae.  Axillae  are  more  easily  visualized  with  increased  system  tilt. 


Reconstructed  results  are  shown  in  Fig.  10.  From  these  preliminary  axillary  imaging  studies,  we  have  shown  that 
axillary  imaging  is  indeed  possible,  even  in  the  case  of  0°  system  tilt  (Fig.  10).  Furthermore,  the  axillae  are  more 
clearly  visualized  when  the  COR  was  offset  from  the  nipple-chest  axis  and  a  modest  amount  of  system  tilt  was  used 
(Fig.  10,  TPB,  60°  and  COMPLEX).  To  simultaneously  image  both  the  breast  and  the  axillae,  with  limited  distortion  in 
the  breast,  an  orbit  such  as  COMPLEX  could  be  employed.  However,  the  COMPLEX  orbit  is  not  practical  because  we 
spent  several  hours — even  with  the  automated  orbit  creation  code  developed — attempting  to  optimize  the  orbit 
considering  both  the  polar  tilt  and  the  ROR  (to  create  COMPLEX,  ~36  hours  were  required  to  customize  the  orbit).  A 
simpler  implementation  method  will  need  to  be  developed  so  that  patients  will  not  be  required  to  lie  on  a  table  for 
extended  periods  of  time. 
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Key  Research  Accomplishments 


Year  2  included  tasks  2(b)-3(a)  from  the  original  Statement  of  Work  (Appendix  A).  Task  2  was  completed  in  its 
entirety.  The  additional  positioning  complication  of  bed  shielding  was  investigated,  preliminary  measurements  of 
system  tilt  for  axillary  imaging  were  completed,  a  program  was  developed  to  post-process  the  acquired  list-mode  data 
to  determine  optimal  acquisition  and  system  parameters,  and  Matlab  programs  were  developed  to  automate  orbit 
creation,  including  dynamic  ROR  control.  Summarizing  specific  progress  for  Task  2: 

■  Experiments  from  the  effects  of  additional  bed  shielding  were  presented  at  the  2 004  IEEE  Nuclear  Science 
Symposium  and  Medical  Imaging  Conference,  16-22  October,  Rome,  Italy,  and  published  in  its  associated 
conference  record  (Appendix  B). 

■  The  results  from  the  preliminary  axillary  imaging  studies  were  submitted  and  accepted  as  an  oral  presentation 
to  the  2 005  IEEE  Nuclear  Science  Symposium  and  Medical  Imaging  Conference,  23-29  October,  San  Juan,  Puerto 
Rico  (Apendix  C).  An  accompanying  manuscript  will  be  submitted  to  the  conference  record. 

*  The  LMPro  and  Matlab  programs  were  developed  and  have  been  included  within  my  dissertation.  An  abstract 
utilizing  results  from  the  LMPro  analysis  were  submitted  and  accepted  to  the  2005  IEEE  Nuclear  Science 
Symposium  and  Medical  Imaging  Conference  (See  Appendix  D  for  submitted  abstract),  in  which  I  am  a  co¬ 
author. 

■  Some  of  the  data  acquired  for  task  3(a)  was  incorporated  into  a  published  article,  "Evaluation  of  Fully  3D 
Emission  Mammotomography  with  a  Compact  Cadmium  Zinc  Telluride  Detector,"  IEEE  Transactions  on 
Medical  Imaging,  2005.  See  attached  article  in  Appendix  E  for  details. 

Related: 

■  Results  from  this  work  were  included  in  my  preliminary  proposal,  "Development  and  Evaluation  of  a 
Dedicated  Emission  Mammotomography  System,"  which  was  accepted  by  my  dissertation  committee  on 
September  16,  2004. 

■  Our  lab  grew  by  three  graduate  students  this  year,  beginning  in  September  2004.  I  have  had  extensive 
opportunities  to  mentor,  train,  and  educate  Priti  Madhav,  Dominic  Crotty,  and  Spencer  Cutler,  for  whom  we 
have  been  co-authors  on  submitted  and  accepted  abstracts  to  the  2005  IEEE  Nuclear  Science  Symposium  and 
Medical  Imaging  Conference.  A  post-doctoral  student.  Dr.  Esther  Dura,  also  joined  the  lab  in  June  2005. 
Already,  we  have  interacted  many  times,  exchanging  ideas  about  image  registration  of  the  functional  and 
structural  reconstructed  images  obtained  from  our  systems. 

■  I  have  been  able  to  apply  knowledge  learned  about  unique  3D  acquisition  orbits  from  this  work  to  X-ray 
computed  mammotomography,  resulting  in  several  published  conference  proceedings  of  the  SPIE.  In  addition,  I 
have  also  contributed  to  a  comparison  of  mammotomography,  which  utilized  the  arbitrary  orbits,  to 
scintimammography. 

■  lam  currently  in  the  process  of  writing  my  dissertation,  which  I  expect  to  defend  within  the  next  few  months. 
All  of  the  results  from  this  grant  have  been  included  within  my  thesis. 

Reportable  Outcomes 

Peer-reviewed  articles 

C.N.  Brzymialkiewicz,  M.P.  Tornai,  R.L.  McKinley,  and  J.E.  Bowsher,  "Evaluation  of  Fully  3D  Emission 

Mammotomography  with  a  Compact  Cadmium  Zinc  Telluride  Detector,"  IEEE  Trans,  on  Med.  Imag.  vol.  24,  no.  7, 

pp.  868-877,  2005. 

M.P.  Tornai,  C.N.  Brzymialkiewicz,  M.L.  Bradshaw,  J.E.  Bowsher,  B.E.  Patt,  J.S.  Iwanczyk,  J.  Li,  and  L.R.  MacDonald, 

"Comparison  of  Compact  Gamma  Cameras  with  1.3mm  and  2.0mm  Quantized  Elements  for  Dedicated  Emission 

Mammotomography,"  IEEE  Trans.  Nucl.  Sci.  (In  Press). 

Presentations,  Published  Abstracts,  and  Published  Proceedings: 
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C.N.  Brzymialkiewicz,  M.P.  Tornai,  R.L.  McKinley,  and  J.E.  Bowsher,  "3D  Data  Acquisition  Sampling  Strategies  for 
Dedicated  Emission  Mammotomography  for  Various  Breast  Sizes,"  Presented  at  the  2004  IEEE  NSS/MIC,  16-22 
October,  Rome,  Italy,  and  published  in  the  2004  IEEE  NSS/MIC  Conference  Record,  vol.  4,  pp.  2596  -  2600. 

M.P.  Tornai,  C.N.  Brzymialkiewicz,  and  S.J.  Cutler  and  P.  Madhav,  "Comparison  of  scintimammography  and 
dedicated  emission  mammotomography,"  Presented  at  the  2004  IEEE  NSS/MIC,  16-22  October,  Rome,  Italy,  and 
published  in  the  2004  IEEE  NSS/MIC  Conference  Record,  vol.  5,  pp.  2818  -  2822. 

M.P.  Tornai,  R.L.  McKinley,  C.N.  Brzymialkiewicz,  P.  Madhav,  S.J.  Cutler,  D.J.  Crotty,  J.E.  Bowsher,  E.  Samei  and 
C.E.  Floyd,  "Design  and  Development  of  a  Fully-3D  Dedicated  X-ray  Computed  Mammotomography  (CmT) 
System,"  Presented  at  the  2005  SPIE  Medical  Imaging  Conference,  San  Diego,  CA,  12  -  17  Feb.  2005,  and  published  in 
Proc.  SPIE:  Physics  of  Medical  Imaging,  5745(1):189-197. 

C.N.  Brzymialkiewicz  and  M.P.  Tornai,  "Novel  Three-dimensional  Acquisition  Orbits  with  a  Dedicated  Emission 
Mammotomography  System,"  Presented  at  the  DOD  Era  of  Hope  Conference,  8-11  June  2005,  Philadelphia,  PA. 

C. N.  Brzymialkiewicz,  R.L.  McKinley,  and  M.P.  Tornai,  "Initial  Patient  Results  from  Dedicated  Emission 

Mammotomography,"  Abstract  submission  accepted  to  the  2005  Nucl.  Sci.  Sump.  &  Med.  Imag.  Conf.,  San  Juan, 
Puerto  Rico,  23-29  Oct.  2005. 

D. J.  Crotty,  C.N.  Brzymialkiewicz,  R.L.  McKinley,  and  M.P.  Tornai,  "Optimizing  Orientation  of  SPECT  and  CT 

Detectors  Through  Quantification  of  Cross  Contamination  in  a  Dual  Modality  Mammotomography  System," 
Abstract  submission  accepted  to  the  2005  Nucl.  Sci.  Symp.  &  Med.  Imag.  Conf.,  San  Juan,  Puerto  Rico,  23-29  Oct. 
2005. 

S.J.  Cutler,  C.N.  Brzymialkiewicz,  and  M.P.  Tornai,  "Investigating  the  Effects  of  Energy  Resolution  in  Dedicated 
Emission  Mammotomography,"  Abstract  submission  accepted  to  the  2005  Nucl.  Sci.  Symp.  &  Med.  Imag.  Conf,  San 
Juan,  Puerto  Rico,  23-29  Oct.  2005. 

P.  Madhav,  C.N.  Brzymialkiewicz,  J.E.  Bowsher,  and  M.P.  Tornai,  "A  Novel  Way  to  Characterize  the  MTF  in  3D  for 
Quantized  SPECT  Cameras  Having  Arbitrary  Trajectories,"  Abstract  submission  accepted  to  the  2005  Nucl.  Sci. 
Symp.  &  Med.  Imag.  Conf,  San  Juan,  Puerto  Rico,  23-29  Oct.  2005. 

Conclusions 

From  these  studies  utilizing  complex  3D  acquisition  orbits,  we  have  demonstrated  that  lesions  near  the  chest  wall  can 
be  visualized  for  both  large  and  small  sized  breasts,  and  further,  that  axillary  imaging  is  possible  with  our  system. 
Acquisition  parameters  were  also  evaluated  to  determine  an  optimal  imaging  protocol  for  potential  future  patient 
studies.  Given  the  very  similar  SNR  and  contrast  values  for  all  examined  orbits,  it  would  seem  that  other  factors  are 
more  important  in  determining  sampling  strategy,  including  the  reconstructed  breast  and  lesion  shape,  visualization  of 
the  area  near  the  chest  wall  (i.e.  lesion  location),  and  area  of  the  reduced-artifact  region  of  the  breast.  However,  one 
could  utilize  a  different  sampling  strategy  depending  on  the  imaging  task.  Future  work  includes  the  addition  of  a 
patient  bed,  which  is  expected  to  complicate  positioning  even  further,  and  which  may  or  may  not  include  lead  lining 
pending  further  studies. 
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APPENDIX  A 


Statement  of  Work 

Task  1  Implement  and  characterize  modified  basis  set  of  orbits  (Months  1-5): 

a.  Implement  dynamic  ROR  capability  in  hardware  to  enable  fully  computer  controlled  gantry.  (Month  1) 

b.  Program  gantry  for:  modified  circle-plus-arc  orbit  with  different  arc  locations,  cloverleaf  orbit.  (Month 

2) 

c.  Perform  phantom  measurements  using  modified  orbits  and  evaluate  projection  data  for  signal-to-noise 
(SNR),  contrast,  and  lesion  visualization  improvement.  (Months  2-5). 

Task  2  Development  and  investigation  of  novel  orbits  (Months  5-22): 

a.  Use  Monte  Carlo  and  Analytic  computer  simulations,  including  computations  of  Orlov  Volumes,  to 
develop  novel  orbits  for  various  breast  shapes  and  sizes  through  parameters  of  viewable  breast  volume, 
radius-of-rotation  (ROR),  camera  tilt  angle  (to  minimize  background  contamination).  (Months  5-12) 

b.  Investigate  the  effect  of  additional  bed  shielding  on  positioning  of  the  compact  gamma  camera. 
(Months  12-15) 

c.  Investigate  acquisition  and  system  parameters,  including  system  tilt  for  axillary  imaging,  total  number 
of  projections  and  hence  angular  sampling,  and  distribution  of  the  scan  time.  (Months  15-18) 

d.  Develop  software  to  automate  process  of  orbit  creation,  utilizing  dynamic  ROR  control.  (Months  18-22) 

Task  3  Experimental  evaluation  of  orbits  (Months  22-36): 

a.  Acquire  experimental  projection  data  using  cold  rod  and  cold  disk  phantoms  for  resolution  and  sampling 
characterization,  respectively.  (Months  22-29) 

b.  Utilize  anthropomorphic  breast  containing  lesions  and  torso  phantoms  to  acquire  projection  data. 
Analyze  reconstructed  images  for  contrast,  signal-to-noise  ratio,  lesion  detectability,  and  extent  of 
artifacts  caused  by  torso  contamination.  (Months  29-36) 
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APPENDIX  B: 

IEEE  2004  Nuclear  Science  Symposium  and  Medical  Imaging  Conference 

Conference  Record 


3D  Data  Acquisition  Sampling  Strategies  for 
Dedicated  Emission  Mammotomography  for  Various 

Breast  Sizes 


Caryl  N.  Brzymialkiewicz,  Member,  IEEE,  Martin  P.  Tomai,  Member,  IEEE, 

Randolph  L.  McKinley,  Member,  IEEE,  Spencer  J.  Cutler,  Member,  IEEE,  and  James  E.  Bowsher,  Member,  IEEE 


Abstract — The  dedicated  emission  mammotomography  system 
developed  in  our  lab  is  in  preparation  for  initial  patient  studies.  As 
a  preliminary  step,  we  evaluate  the  effect  of  breast  size  and  lesion 
location  on  this  paradigm.  The  hemispherical  positioning  gantry 
allows  ample  flexibility  in  sampling  a  pendant,  uncompressed 
breast.  Recently  acquired,  realistic  anthropomorphic  torso  (which 
includes  the  upper  portion  of  the  arm)  and  breast  phantoms 
emphasize  the  necessity  of  employing  unique  camera  trajectories 
(orbits)  rather  than  simple  VAOR  camera  trajectories.  Several 
novel  3D  orbits  have  been  implemented  with  fully  contoured 
radius-of-rotation  capability  to  compensate  for  the  positioning 
demands  that  are  required  for  different  breast  sizes.  While  a 
general  orbit  design  may  remain  the  same  between  two  different 
breasts,  the  absolute  polar  tilt  range  and  ROR  range  may  vary. 
We  have  demonstrated  that  with  increased  polar  camera  tilt, 
employing  3D  data  acquisition  camera  trajectories,  lesions  near 
the  chest  wall  can  be  visualized  for  both  large  and  small  sized 
breasts. 

Index  Terms — Breast  imaging,  emission  mammotomography, 
single  photon  emission  computed  tomography,  three-dimensional 
orbits,  CZT  detector 

I.  Introduction 

Our  group  has  previously  investigated  several  parallel-beam 
fully  3D  acquisition  orbits  about  a  pendant  breast,  both  with 
and  without  additional  torso  backgrounds.  Results  thus  far  have 
shown  that  no  single  camera  trajectory  will  equally  optimize 
image  quality  for  various  lesion  sizes  located  in  different  areas 
of  the  breast  [1] — [3],  However,  some  orbits  may  be  more  easily 
adaptable  to  varieties  of  breast  shapes  and  sizes.  For  this  set 
of  experiments,  we  focus  on  three  specific  orbit  types:  Vertical 
Axis  of  Rotation  (VAOR),  Tilted  Parallel  Beam  (TPB),  and 
Projected  Sine  Wave  (PROJSINE). 

The  dedicated  emission  mammotomography  system  we  have 
been  developing  in  our  lab  is  in  preparation  for  initial  patient 
studies.  As  a  preliminary  step,  we  evaluate  the  effect  of 
various  breast  sizes  and  lesion  location  on  this  paradigm.  The 
hemispherical  positioning  gantry  allows  ample  flexibility  in 
sampling  a  pendant,  uncompressed  breast.  The  gantry  itself 

Manuscript  submitted  November  1,  2004.  This  work  was  supported  by  the 
National  Institutes  of  Health  under  grant  R01-CA96821  and  by  U.S.  Army 
Grant  No.  DAMD17-03-1-0558. 

C.N.  Brzymialkiewicz,  M.P.  Tomai,  R.L.  McKinley,  and  S.J.  Cutler  are  with 
the  Departments  of  Biomedical  Engineering  and  Radiology;  J.E.  Bowsher  is 
with  the  Department  of  Radiology,  Duke  University  Medical  Center,  Durham, 
NC  27710  USA  (email:cnb3@duke.edu). 


has  been  described  in  detail  elsewhere  [1],  [2].  The  compact 
gamma  camera  currently  attached  to  the  system  is  the  CZT- 
based  Lumagem  3200-S  ( Gamma  Medica,  Northridge,  CA) 
with  2.5  mm2  discrete  pixels,  with  a  measured  energy  resolution 
of  6.1%  FWHM  at  140  keV  and  sensitivity  of  37.9  cps/MBq. 

Key  benefits  of  utilizing  3D  acquisition  orbits  include  the 
ability  to  contour  the  breast,  reducing  distance-dependent  spa¬ 
tial  resolution  limitations,  and  the  capability  of  imaging  further 
into  the  breast,  near  the  chest  wall,  by  employing  camera  polar 
tilt.  Recently  acquired,  realistic  anthropomorphic  torso  (which 
includes  the  upper  portion  of  the  arm)  and  breast  phantoms 
emphasize  the  necessity  of  employing  unique  camera  trajecto¬ 
ries  (orbits)  rather  than  simple  VAOR  camera  trajectories.  For  a 
simple  circular  orbit  (VAOR,  for  example,  with  360°  azimuthal 
rotation,  no  polar  camera  tilt),  the  amount  of  viewable  breast 
volume  is  limited  by  the  placement  of  the  camera  next  to 
the  torso.  Due  to  physical  limitations,  the  tissue  immediately 
anterior  to  the  chest  wall  is  not  within  the  field-of-view  (FOV) 
(Fig.  1,  VAOR,  white  arrow).  By  simply  increasing  the  camera 
polar  tilt  angle  to  45°,  the  center-of-rotation  (COR)  moves 
further  into  the  breast  as  compared  to  VAOR,  thus  allowing 
the  entire  volume  to  be  viewed.  However,  artifacts  result  from 
the  incompletely  sampled  TPB  orbit.  Hence,  more  complex 
3D  camera  trajectories  which  more-nearly-completely  sample 
the  breast  while  minimizing  direct  lines-of-sight  of  the  heart 
and  liver  are  desired.  Furthermore,  complex  3D  orbits  can  be 
designed  to  avoid  physical  limitations  imposed  by  the  patient 
(e.g.  arms)  or  bed. 

One  type  of  complex  3D  orbit  investigated  in  this  work  is 
the  PROJSINE  (Fig.  2),  which  is  a  sine  wave  projected  onto 
the  hemisphere,  where  the  number  of  ‘lobes’  and  the  polar  tilt 
range  can  be  varied.  Three  lobes  minimize  the  number  of  times 
the  camera  must  reach  <p= 0°  (a  position  that  requires  larger 
radius-of-rotation  (ROR)).  Two  variations  of  PROJSINE  were 
examined  within  this  work,  one  going  from  0°  to  45°,  and  the 
other  from  15°  to  45°.  By  allowing  the  camera  to  start  at  a 
polar  tilt  angle  of  15°  rather  than  0°,  the  COR  moves  further 
into  the  breast,  yielding  an  increased  viewable  breast  volume 
including  imaging  near  the  chest  wall,  and  decreased  maximum 
ROR. 
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Fig.  1.  COR  moves  further  into  the  breast  with  increased  camera  polar  tilt, 
allowing  for  imaging  of  the  breast  near  the  chest  wall. 


Fig.  3.  Photographs  of  the  1 060  ml  (left)  and  325  ml  (right)  anthropomorphic 
breast  phantoms.  Dimensions  of  the  breast  phantoms  (given  in  [cm]  as: 
superior-inferior  x  nipple-chest  x  medial-lateral)  for  the  1060  ml  phantom: 
15.9  x  12.1  x  17.8;  for  the  325  ml  phantom:  12.9  x  4.4  x  13.5.  Two  sets  of 
embedded  lesions  were  inserted  into  each  phantom,  with  lesion  sizes  as  labeled. 


TABLE  I 


Fig.  2.  Within  each  box,  the  left  plot  shows  polar  camera  tilt  as  a  function 
of  the  azimuthal  angle;  the  right  plot  shows  the  radius  of  rotation,  shown  here 
for  the  1060  ml  breast  phantom. 


Parameters  Used  for  each  Breast-Only  Acquisition,  with  all 
ORBITS  ACQUIRED  OVER  A  360°  AZIMUTHAL  RANGE  (0). 


Orbit 

(Range, 

min:max) 

32; 

COR 
Location 
(relative 
to  the 
nipple)* 
(cm) 

i  ml 

ROR* 

(Range, 

mimmax) 

(cm) 

106 

COR 
Location 
(relative 
to  the 
nipple)* 
(cm) 

0  ml 

ROR* 

(Range, 

mimmax) 

(cm) 

VAOR 

0° 

-5.5 

6.1 

-1.0 

8.0 

TPB 

45° 

-2.6 

1.1 

1.7 

4.0 

PROJSINE 

MPlfT  flilP 

-5.5 

0.0  :  7.0 

-1.0 

0.8  :  9.4 

PROJSINE 

■  LWKM 

-2.6 

1.7 

3.1  :  7.7 

1 Negative  values  indicate  CORs  outside  and  below  the  nipple. 
t  ROR  defined  as  the  distance  from  the  center  of  the  detector  to  the  COR. 


II.  Methods 

A.  Isolated  Breast  Phantom 

Two  custom  breast  phantoms  (Fig.  3)  ( Radiology  Support 
Devices,  Newport  Beach,  CA)  were  used  in  conjunction  with 
the  3D  orbits.  Three  lesions  were  placed  in  each  breast  (Fig. 
3);  one  0.45  ml  catheter  balloon  lesion  ( Harvard  Apparatus, 
Holliston,  MA)  was  glued  onto  the  anterior  chest  wall,  one  0.43 
ml  acrylic  sphere  lesion  ( Data  Spectrum  Corp.,  Hillsborough, 
NC)  was  placed  as  centrally  in  the  breast  as  possible,  and 
one  0.26  ml  lesion  ( Data  Spectrum  Corp.)  was  placed  laterally 
within  the  breast.  The  absolute  activity  for  the  lesions  in  the 
325  ml  breast  was  21.1  /iCi/ml;  in  the  1060  ml  breast,  the 
lesion  activity  was  20.0  /rCi/ml.  Ten  minute  scans  were  per¬ 
formed,  with  acquisition  times  increased  in  subsequent  scans  to 
account  for  radioactive  decay.  A  ±8%  energy  window  about  the 
140  keV  photopeak  of  Tc-99m  was  used.  Three  lesiombreast 
concentration  ratios  were  used,  including  10:1,  6:1,  and  3:1. 
Specific  orbit  parameters  used  for  each  scan  are  given  in  Table 
I.  VAOR  and  TPB  each  had  128  projections;  the  PROJSINE 
orbits  had  150  projections. 
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V 

(Range, 

minrmax) 


COR  Location 
(relative  to 
the  nipple)* 
(cm) 


-1.0 


1.7 


-1.0 


1.7 


Negative  values  indicate  CORs  outside  and  below  the  nipple. 

* ROR  defined  as  the  distance  from  the  center  of  the  detector  to  the  COR. 


VAOR 


TPB 


PROJSINE 


PROJSINE 


B.  Breast  Phantom  with  Filled  Torso  Background 

Two  lesions  were  embedded  in  the  1060  ml  breast  phantom, 
which  was  then  attached  to  a  filled  anthropomorphic  torso  phan¬ 
tom  that  included  a  heart,  lungs,  and  liver  ( RSD )  (Fig.  4).  A  1 .07 
ml  (M  .27  cm  diam.)  lesion  was  placed  near  the  chest  wall;  a 
2.27  ml  (~1.6  cm  diam.)  lesion  (RSD)  was  placed  in  the  center 
of  the  breast.  The  absolute  lesion  activity  was  10.4  /rCi/ml;  the 
lesion:torso-and-breast:heart-and-liver  concentration  ratio  was 
11:1:10.  Projection  data  were  acquired  for  20  minutes,  with 
subsequent  scans  adjusted  to  compensate  for  radioactive  decay. 
A  ±8%  energy  window  was  again  used.  Scans  were  acquired 
with  and  without  the  1/4”  lead  sheet,  which  was  securely  taped 
to  the  phantom  to  simulate  shielding  of  the  ‘patient’  along  with 
the  fully  pendant  exposed  breast  of  interest. 


0.43  ml  0.26  ml 


Fig.  5.  Top  three  rows  are  images  obtained  from  reconstructed  data  (OSEM, 
2.5mm  voxels,  2nd  iteration  shown,  attenuation  correction,  3  summed  slices 
to  reduce  noise)  for  the  1060  ml  breast  phantom,  with  a  10:1  concentration 
ratio,  and  showing  all  lesions.  Note  that  the  0.45  ml  lesion  (located  near  the 
chest  wall)  is  not  visible  with  VAOR.  Bottom  row  is  the  Maximum  Intensity 
Projection  (MIP)  through  the  breast  volume. 


for  the  breast  only  studies;  a  9-pixel  ROI  was  drawn  inside  and 
about  the  2.27  ml  lesion  for  the  breast+torso  studies. 

The  SNR  was  measured  as  the  difference  in  mean  pixel 
values  of  the  lesion  and  breast  background,  divided  by  the 
standard  deviation  of  the  uniform  background  ((ROI;es,on- 
ROI,„.,„ )l  t7  ijr,j:,.f) .  where  ROFesjOTl  and  ROI  breast  4**®  the 
mean  values  of  the  counts  within  each  respective  re¬ 
gion  and  abreast  is  the  standard  deviation  of  the  uniform 
breast  background.  Lesion  contrast  is  defined  as  the  sig¬ 
nal  difference  divided  by  the  background  signal  ((ROI(esjora- 
ROI(,reast)/ROI(,reasf). 


C.  Image  Reconstruction  and  ROI  Analysis 

Images  were  reconstructed  from  projection  data  using  an 
OSEM  algorithm  with  ray-tracing  capability,  with  8  subsets, 
5  iterations,  and  calculated  attenuation  correction  assuming  a 
uniform  emission  volume  [4],  [5].  For  ease  of  display  only 
(because  the  lesions’  locations  are  through  various  planes  of  the 
breast).  Maximum  Intensity  Projections  (MIP)  are  used  to  sum¬ 
marize  the  reconstructed  data.  Regions-of-interest  (ROIs)  were 
drawn  in  the  reconstructed  images  both  inside  the  lesions  and 
circumferentially  about  each  lesion  in  the  breast  background. 
A  9-pixel  ROI  was  drawn  inside  and  about  the  0.43  ml  lesion 


III.  Results  &  Discussion 
As  seen  in  Fig.  5,  the  lesion  nearest  the  chest  wall  (0.45  ml) 
cannot  be  seen  with  VAOR.  This  lesion  is  not  visible  for  the 
1060  ml  breast  nor  the  325  ml  breast  with  VAOR  (Figs.  6  and 
7),  a  result  of  the  physical  limitation  of  placing  the  detector  near 
the  torso.  For  the  3:1  concentration  ratio,  none  of  the  lesions 
were  visible  in  the  1060  ml  breast  (Fig.  6),  while  all  lesions 
were  visible  in  the  325  ml  breast  phantom  (Fig.  7)  (with  the 
exception  of  the  chest  wall  lesion  using  VAOR).  This  may  be 
expected  given  the  increased  ROR  with  the  larger  breast  and 
increased  amount  of  attenuating  tissue. 
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Fig.  7.  Maximum  Intensity  Projections  (Sagittal  view)  for  the  325  ml  breast 
for  all  concentration  ratios  and  orbits,  as  labeled.  White  arrow  points  to  cold 
area,  due  to  nylon  screw. 


SNR 


Fig.  8.  SNR  and  Contrast  results  for  the  centrally-located  0.43  ml  lesion  (2nd 
iteration)  embedded  in  the  1060  ml  breast. 


Fig.  10.  Images  from  reconstructed  data  (OSEM,  2.5mm  voxels,  2nd  iteration 
shown,  attenuation  correction,  3  summed  slices,  transaxial  view)  for  the  1060 
ml  breast  and  filled  torso,  with  and  without  the  lead  sheet.  Similar  to  the 
breast-only  studies,  the  smaller  lesion  is  not  visible  with  VAOR  and  artifacts 
move  further  into  the  breast  for  the  PROJSINE  15-45°  case  as  compared  to 
PROJSINE  0-45°. 


As  expected,  the  TPB  orbit  yields  a  distorted  breast  shape 
due  to  incomplete  sampling.  Indeed,  the  lesions  are  also  some¬ 
what  distorted  with  TPB.  The  more-nearly-completely  sampled 
PROJSINE  orbits  recovers  the  breast  shape  more  accurately 
than  TPB,  however,  artifacts  near  the  chest  wall  are  still  visible. 
Note  that  there  is  no  activity  “above”  the  breast  volume  (i.e.  in 
the  chest  wall  region),  yet  because  of  the  sampling,  in  which  all 
activity  is  not  within  the  FOV  for  all  vantages  of  the  camera, 
the  reconstructed  images  have  incorrectly  placed  activity.  By 
sacrificing  some  polar  sampling  of  PROJSINE  by  starting  at 
15°  polar  tilt  rather  than  0°,  the  COR  moves  further  into  the 
breast,  thus  increasing  the  reduced-artifact  region  of  the  breast 
volume  and  allowing  easier  visualization  of  the  lesion  near  the 
chest  wall  (Fig.  6,  white  arrows). 


0-7803-8700-7/04/$20.00  (C)  2004  IEEE 


2599 


Without  lead  sheet  With  lead  sheet 


IV.  Conclusions 


on  •'t 


Fig.  11.  Images  from  reconstructed  data  (OSEM,  2.5  mm  voxels,  2nd  iteration 
shown,  attenuation  correction,  summed  slices,  sagittal  view)  for  the  1060  ml 
breast  and  filled  torso,  with  and  without  the  lead  sheet.  More  background 
contamination  is  visible  without  the  lead  sheet;  increased  tilt  views  the  heart 
and  liver  more  directly. 


SNR 


**#  '  No  shield,  VAOR 

♦  No  shield.  TPB  45 

♦  No  shield,  Projsiiic  0-45 

'  No  shield,  Projsine  15-45 
^  With  shield,  VAOR 
HQ  With  shield,  1PB  45 
'0  With  shield,  Projsine  0-45 
With  shield,  Projsine  15-45 

Fig.  12.  SNR  and  Contrast  results  for  the  2.27  ml  lesion.  Iterations  labeled 
next  to  TPB  45°.  Patient  shielding  yields  slightly  lower  SNR  and  contrast 
values. 


For  each  concentration  ratio,  all  orbits  yield  similar  SNR  and 
contrast  values  for  the  1060  ml  breast  phantom  (Fig.  8);  for  the 
325  ml  breast,  TPB  45°  seems  to  outperform  the  others  (Fig. 
9). 

Reconstructed  results  for  the  shielding  experiments  are 
shown  in  Figs.  10  and  11.  No  patient  shielding  yields  slightly 
higher  SNR  and  contrast  values  than  with  patient  shielding 
(Fig.  12).  We  intend  to  repeat  the  study  with  breast-only  filled 
measurements,  prior  to  attaching  the  breast  to  the  filled  torso 
phantom. 


Orbits  that  include  increased  camera  polar  tilt  angles  allow 
for  imaging  of  lesions  near  the  chest  wall  and/or  smaller  breasts. 
Because  previous  studies  have  suggested  the  tomographic  sam¬ 
pling  criteria  for  completeness  may  not  be  as  strict  as  initially 
anticipated,  we  relaxed  the  restrictions  on  the  camera  polar  tilt. 
This  allowed  the  COR  to  move  further  into  the  breast,  which 
increased  the  reduced-artifact  region  of  the  breast,  thus  allowing 
easier  visualization  of  lesions  near  the  chest  wall.  These  studies 
were  intentionally  not  performed  with  clinical  count  densities 
since  we  wanted  to  determine  the  effect  of  the  sampling 
strategy  itself.  Given  the  very  similar  SNR  and  contrast  values 
for  all  examined  orbits,  it  would  seem  that  other  factors  are 
more  important  in  determining  sampling  strategy,  including 
the  reconstructed  breast  and  lesion  shape,  visualization  of  the 
area  near  the  chest  wall  (i.e.  lesion  location),  and  area  of  the 
reduced-artifact  region  of  the  breast.  However,  one  could  utilize 
a  different  sampling  strategy  depending  on  the  imaging  task. 
For  example,  if  a  lesion  is  suspected  near  the  chest  wall,  then 
to  confirm  diagnosis  of  a  lesion  in  that  region  one  may  utilize 
more  camera  tilt;  to  determine  whether  a  smaller  lesion  and/or 
lesion  of  low  radiopharmaceutical  uptake  is  present  anywhere 
in  the  breast,  an  orbit  such  as  TPB,  with  its  decreased  ROR, 
may  be  more  appropriate.  Further  studies  will  be  necessary  to 
fully  characterize  the  sensitivity  of  the  system  in  various  breast 
sizes. 
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Initial  Patient  Results  from  Dedicated  Emission  Mammotomography 
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A  patient  bed  was  recently  developed  and  subsequent  modifications  made  for  3D,  pendant  breast  imaging  with  a 
dedicated  emission  mammotomography  system  that  incorporates  a  discretized  16x20cmA2  CZT  gamma  camera  with 
2.5mm  pixels.  Goals  for  these  experiments  include:  (1)  designing  acquisition  orbits  to  minimize  effects  from  complications 
that  arise  with  patient  positioning;  (2)  determining  the  optimal  amount  of  system  tilt  to  adequately  image  the  breast;  (3) 
determining  whether  axillary  imaging  is  possible  in  this  configuration.  Data  was  initially  acquired  with  Tc-99m  filled 
anthropomorphic  breast  and  torso  phantoms  for  each  of  three  cases:  activity  in  the  breast  only;  activity  in  the  breast  and 
torso;  activity  in  both,  with  lead  shielding  placed  on  the  pallet  (to  mimic  an  intended  custom-designed  lead-lined  pallet). 
Projection  data  were  acquired  for  20  minutes  using  both  simple  and  complex  3D  orbits.  Optimal  positioning  parameters 
determined  from  the  phantom  studies  are  then  applied  to  a  patient  scan.  We  have  found  that  acquisition  orbits  can  be 
modified  to  accommodate  the  bed's  physical  hindrances  and  still  recover  most  of  the  volume  of  interest,  while 
minimizing  distortion  and  artifacts  in  the  reconstructed  images.  Initial  patient  results  with  the  mammotomography 
system  are  presented.  Overall,  the  results  indicate  this  approach  could  be  a  beneficial  diagnostic  tool  to  aid  in  the 
characterization  of  lesions. 

This  work  was  funded  by  DOD  DAMD17-03-1-0558,  NIH  R01-CA96821,  and  in  part  by  W81XWH-05-1-0280. 
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Investigating  the  Effects  of  Energy  Resolution  in 
Dedicated  Emission  Mammotomography 

Spencer  J.  Cutler1,2,  Caryl  N.  Brzymialkiewicz1,2,  Martin  P.  Tornai1,2 
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This  study  probes  the  recent  debate  over  the  necessity  for  good  energy  resolution  for  breast  lesion  imaging  with 
dedicated  single  photon  emission  mammotomography,  an  uncompressed  breast,  fully-3D  imaging  modality  developed  in 
our  lab.  The  system  consists  of  a  commercial,  discretized  CZT  gamma  camera  having  6.7%  FWHM  energy  resolution  (at 
140keV)  and  intrinsic  spatial  resolution  corresponding  to  the  2.5mm  square  pixel,  and  is  used  on  a  fully-3D  positioning 
gantry.  The  energy  resolution  on  this  high  performance  imaging  device  can  be  effectively  degraded  by  using  wider 
energy  windows  in  list  mode  acquired  data.  Multiple  3D  orbits  about  Tc-99m  filled  anthropomorphic 
breast  and  torso  phantoms  are  acquired.  The  full  spectrum  list  mode  data  files  were  multiply  processed  to  obtain  images 
of  varying  energy  window  widths  (from  symmetric  6%  to  an  asymmetric  18%  (-12+6))  but  with  the  same  projection  image 
count  density.  Counts  were  randomly  subsampled  from  the  entire  list  mode  data  set  in  order  to  maintain  an  equivalent 
level  of  noise.  The  data  was  then  reconstructed  using  OSEM  for  various  iterations.  Regions  of  interest  were  drawn  in  and 
about  spherical  lesions  of  different  volumes  to  determine  signal-to-noise  ratios  and  contrasts  for  each  iteration.  Results 
clearly  illustrate  both  visual  and  quantitative  differences  between  the  various  energy  windows,  with  smaller  energy 
windows  (better  ERes)  having  better  image  quality.  Results  for  clinical  count  densities  are  also  investigated. 

This  work  was  funded  by  NIH  R01-CA96821,  in  part  by  DOD  DAMD17-03-1-0558,  and  Duke  BME. 
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Evaluation  of  Fully  3-D  Emission 
Mammotomography  With  a  Compact 
Cadmium  Zinc  Telluride  Detector 
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Abstract — A  compact,  dedicated  cadmium  zinc  telluride  (CZT) 
gamma  camera  coupled  with  a  fully  three-dimensional  (3-D) 
acquisition  system  may  serve  as  a  secondary  diagnostic  tool  for 
volumetric  molecular  imaging  of  breast  cancers,  particularly  in 
cases  when  mammographic  findings  are  inconclusive.  The  devel¬ 
oped  emission  mammotomography  system  comprises  a  medium 
field-of-view,  quantized  CZT  detector  and  3-D  positioning  gantry. 
The  intrinsic  energy  resolution,  sensitivity  and  spatial  resolution 
of  the  detector  are  evaluated  with  Tc-99m  (140  keV)  filled  flood 
sources,  capillary  line  sources,  and  a  3-D  frequency-resolution 
phantom.  To  mimic  realistic  human  pendant,  uncompressed 
breast  imaging,  two  different  phantom  shapes  of  an  average  sized 
breast,  and  three  different  lesion  diameters  are  imaged  to  evaluate 
the  system  for  3-D  mammotomography.  Acquisition  orbits  not 
possible  with  conventional  emission,  or  transmission,  systems  are 
designed  to  optimize  the  viewable  breast  volume  while  improving 
sampling  of  the  breast  and  anterior  chest  wall.  Complications  in 
camera  positioning  about  the  patient  necessitate  a  compromise  in 
these  two  orbit  design  criteria.  Image  quality  is  evaluated  with 
signal-to-noise  ratios  and  contrasts  of  the  lesions,  both  with  and 
without  additional  torso  phantom  background.  Reconstructed 
results  indicate  that  3-D  mammotomography,  incorporating  a 
compact  CZT  detector,  is  a  promising,  dedicated  breast  imaging 
technique  for  visualization  of  tumors  <  1  cm  in  diameter.  Addition¬ 
ally,  there  are  no  outstanding  trajectories  that  consistently  yield 
optimized  quantitative  lesion  imaging  parameters.  Qualitatively, 
imaging  breasts  with  realistic  torso  backgrounds  (out-of-field 
activity)  substantially  alters  image  characteristics  and  breast 
morphology  unless  orbits  which  improve  sampling  are  utilized.  In 
practice,  the  sampling  requirement  may  be  less  strict  than  initially 
anticipated. 

Index  Terms — Breast  imaging,  CZT  detector,  emission  imaging, 
mammotomography,  single  photon  emission  computed  tomog¬ 
raphy,  SPECT,  three-dimensional  orbits. 
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I.  Introduction 

COMPACT,  dedicated  gamma  cameras  coupled  with  a 
fully  three-dimensional  (3-D)  acquisition  system  can 
serve  as  a  secondary  diagnostic  tool  for  imaging  of  breast 
cancers,  particularly  in  cases  when  mammographic  findings 
are  inconclusive.  The  significance  of  a  compact  camera  in 
single  photon  emission  imaging  is  its  ability  to  facilitate  closer 
proximity  imaging,  thus  minimizing  distance-dependent  spatial 
resolution  limitations.  An  advantage  of  employing  nuclear 
medicine  techniques  stems  from  the  fact  that  the  accuracy  of 
scintimammography  is  unaffected  by  breast  density:  studies 
have  shown  women  with  dense  breasts  may  benefit  most  from 
this  complementary  imaging  technique  [1],  [2].  Thus,  several 
groups  have  investigated  compact  position-sensitive  photo¬ 
multiplier  tube  (PSPMT)-based  gamma  cameras  for  the  breast 
imaging  paradigm  [3]— [6]. 

Compared  to  planar  imaging  alone,  there  is  a  significant 
increase  in  the  accuracy  of  breast  imaging  for  the  diagnosis 
of  primary  breast  cancer  when  mammography  and  scintimam¬ 
mography  are  combined  [7].  Diagnostic  accuracy  may  be 
further  improved  when  using  tomography  rather  than  a  simple 
combination  of  scintimammography  and  mammography,  given 
that  contrast  improves  with  such  a  3-D  technique.  In  addition, 
the  3-D  localization  of  the  tumor  is  provided,  as  is  the  potential 
for  in  vivo  quantification  of  lesion  activity.  Single  photon 
emission  computed  tomography  (SPECT)  specifically  for 
breast  imaging  has  been  investigated  both  with  clinical  gamma 
cameras  [8]— [17]  and  dedicated  systems  [4],  [17]— [22].  Our 
group  has  previously  demonstrated  that  a  compact,  fully  3-D 
hemispherical  positioning  gantry  that  allows  for  imaging  with 
a  detector  positioned  anywhere  about  a  pendant,  uncompressed 
breast  phantom  significantly  improves  image  quality  and  lesion 
visualization  [20],  [21],  Those  studies  were  performed  with 
quantized  Nal(Tl)  scintillator,  PSPMT-based  compact  gamma 
cameras. 

While  these  dedicated  PSPMT  cameras  offer  an  improvement 
over  clinical  gamma  cameras,  the  energy  discrimination  can  be 
further  improved  by  using  cadmium  zinc  telluride  (CZT)  semi¬ 
conductor  detectors.  The  use  of  CZT  detectors  for  the  detection 
of  breast  lesions  has  been  previously  proposed  and  investigated 
through  SPECT  simulation  studies  [22]  and  preprototype  mea¬ 
surements  [18],  and  scintimammography  measurements  [23]. 
Solid-state  detectors  made  of  CZT  boast  a  markedly  improved 
energy  resolution  over  scintillator-based  Nal(Tl),  with  reported 
measurements  of  4%-7%  full- width  at  half-maximum  (FWHM) 
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at  140.6  keV  [23]— [27],  as  compared  to  NaI(Tl)’s  typical  en¬ 
ergy  resolution  of  9%-10%  [28].  Indeed,  for  compact  gamma 
cameras  utilizing  quantized  elements,  the  energy  resolution,  and 
hence  scatter  rejection,  is  usually  poorer  still.  Scintimammog- 
raphy  phantom  imaging  studies  indicate  that  scatter  from  car¬ 
diac  and  hepatic  background  can  significantly  contaminate  the 
primary  lesion  signal  [29].  While  it  is  generally  accepted  that 
better  energy  resolution  performance  yields  improved  emission 
image  quality,  it  is  unclear  what  specific  effect  this  will  have  in 
the  3-D  breast  imaging  paradigm  [20],  [21], 

Using  a  fully  3-D  acquisition  geometry  discussed  here,  there 
is  an  infinite  set  of  camera  trajectories,  or  orbits,  with  which 
to  acquire  projection  data.  The  camera  motions  and  techniques 
presented  here  can  also  directly  be  applied  to  other  imaging 
modalities,  such  as  cone-beam  X-ray  transmission  imaging. 
The  gantry’s  several  degrees  of  freedom  allow  any  detector, 
source/detector  or  even  detector/detector  connected  to  it  to 
be  positioned  anywhere  about  a  2-7T  steradian  solid  angle  of 
coverage.  For  example,  by  replacing  the  gamma  camera  utilized 
here  with  a  flat-panel  digital  X-ray  detector,  combined  with  an 
opposing  X-ray  tube,  cone-beam  transmission  acquisitions  of  a 
pendant  breast  phantom  can  be  performed  [30],  [31]. 

While  a  simple,  vertical-axis-of-rotation  orbit  for  an  emis¬ 
sion  imaging  camera  [13]  with  a  parallel-hole  collimator 
yields  a  completely  sampled  volume  [32],  [33],. the  imaged 
breast  volume,  for  example,  is  smaller  than  for  an  acquisition 
orbit  which  includes  camera  tilt.  Due  to  positioning  restraints 
imposed  by  the  detector,  the  furthest  into  the  breast  one  can 
image  without  tilt  is  limited  to  the  plane  perpendicular  to  the 
nipple-chest  axis  where  the  camera  is  adjacent  to  the  torso. 
Cameras  with  a  large  dead  edge  further  decrease  the  viewable 
breast  volume.  One  main  advantage  of  fully  3-D  orbits  that  uti¬ 
lize  increased  polar  tilt  is  an  increased  viewable  breast  volume, 
potentially  into  the  chest  wall  region  [20].  However,  with 
increased  tilt  angle,  the  camera  may  directly  view  cardiac  and 
hepatic  regions  that  take  up  the  injected  radiopharmaceutical, 
which  could  contribute  primary  background  contamination  to 
the  2-D  projection  images  of  the  breast. 

Using  a  compact  CZT  detector  to  image  lesions  in  a  uni¬ 
form  breast  background,  the  primary  aim  of  this  study  is 
to  evaluate  the  imaging  ability  of  the  developed  emission 
mammotomography  system,  both  with  and  without  additional 
torso  backgrounds,  with  some  previously  investigated  and 
newer  trajectories.  To  our  knowledge,  this  is  the  first  imple¬ 
mentation  and  evaluation  of  a  fully  developed  CZT  detector 
and  system  for  dedicated  tomographic  imaging  of  the  breast 
(mammotomography) . 

H.  Planar  and  SPECT  Performance  Characteristics 
of  the  CZT  Detector 

A.  Methods 

The  commercially  available,  production  LumaGEM™ 
3200-S  camera  ( Gamma  Medica  Inc.,  Northridge,  CA)  (Fig.  1) 
utilizes  a  60  x  84  array  of  2.5  x  2.5  x  6  mm3  quantized  CZT 
elements.  The  detector  is  kept  at  15  °C,  using  water  cooling. 
We  evaluated  the  camera  with  an  unmatched  parallel  beam, 
lead  collimator  with  hexagonally  arranged  holes  (1 .22-mm  hole 
size  flat-to-flat,  0.2-mm  septa,  25.4-mm  height). 


Fig.  1.  Photograph  of  the  compact  CZT  detector  with  16  x  20  cm2  field  of 
view  (FOV). 


With  the  collimator  removed,  the  intrinsic  energy  resolution 
was  measured  for  10  M  total  counts  using  a  point  source  of 
Tc-99m  located  approximately  1  m  from  the  camera  surface. 
Count  rates  were  ~17.5  kcps  using  an  open  energy  window 
(~10  kcps  in  a  16%  symmetric  (±8%)  window  about  the  pho¬ 
topeak).  Sensitivity,  defined  as  (counts/sec)/mCi,  was  measured 
with  a  14-cm-diameter  petri  dish  filled  with  an  ~3-mm  layer 
of  2.6  mCi  of  Tc-99m,  placed  10  cm  from  the  collimator  face. 
Counts  were  recorded  for  60  s.  The  dose  calibrator  in  our  lab 
has  a  ±5%  error.  The  sensitivity  measurement  is  based  on  a 
16%  symmetric  (±8%)  energy  window  about  the  140-keV 
photopeak,  which  is  the  window  size  used  for  all  of  the  studies 
performed  to  characterize  the  detector  (all  those  found  in 
Section  H).  Energy  windows  are  applied  globally  to  all  pixels, 
i.e.,  there  is  not  a  separate  energy  window  for  each  discrete 
pixel,  though  in  principle  there  could  be. 

Planar  spatial  resolution  was  measured  by  acquiring  images 
of  a  single  capillary  tube  in  air,  filled  with  3.8  mCi  of  Tc-99m  in 
6.0  cm  of  its  length.  Projection  data  were  acquired  with  the  line 
source  moved  across  the  detector  face  both  horizontally  and  ver¬ 
tically  with  respect  to  the  camera,  in  1-mm  increments  for  a  total 
of  10  mm.  For  both  the  horizontal  and  vertical  directions,  the 
tube  was  then  placed  at  the  position  yielding  the  lowest  FWHM, 
and  projection  data  were  acquired  at  a  collimator-to-source  dis¬ 
tance  of  1  through  10  cm.  The  expected  system  spatial  resolu¬ 
tion  is  a  function  of  the  intrinsic  and  collimator  resolutions.  Col¬ 
limator  resolution  (i?coii)  was  calculated  by  i?coii  —  [d(loff  + 
b  +  c)]//eflF,  where  d  is  the  hole  diameter,  lcft  is  the  ‘effective 
length’  of  the  collimator  holes,  b  is  the  distance  from  the  colli¬ 
mator  to  the  radiation  source,  and  c  is  the  collimator-to-detector 
distance  [34],  [35].  The  effective  collimator  length  is  calculated 
by  Zeff  =  l  —  2 /z-1,  where  l  is  the  length  of  the  holes  and  /j 
is  the  linear  attenuation  coefficient  of  the  collimator  material 
(lead)  [34], 

Reconstructed  spatial  resolution  was  determined  from 
SPECT  measurements  of  two  capillary  tubes  in  air,  one  filled 
with  3.8  mCi  of  Tc-99m  in  6.0  cm  of  its  length,  the  other  filled 
with  1.8  mCi  of  Tc-99m  in  3.0  cm,  both  suspended  along  the 
vertical  axis  of  rotation  with  one  placed  on  the  axis  of  rotation 
and  one  offset  by  2  cm.  Simple  circular  tomography  over 
360°  was  used  for  these  measurements  with  128  projections, 
at  5  s/projection.  Acquisitions  were  made  at  3-,  5-,  and  7-cm 
rotation  radii.  Note  that  the  radius-of-rotation  (ROR)  is  defined 
as  the  perpendicular  distance  from  the  center  of  the  detector 
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Fig.  2.  (Left)  Photograph  of  the  mini-rod  phantom,  with  diameters  labeled 
(in  mm),  with  a  pitch  equal  to  twice  their  diameter.  (Middle)  OSEM 
reconstructed  mini-rod  data  (8  subsets,  20  iterations,  2.5-mm  voxels,  5  summed 
slices  to  reduce  noise,  correction  for  COR  partial  pixel  shift)  acquired  with 
4.4-cm  ROR.  (Right)  Profile  (one-pixel  wide)  drawn  through  the  3.1 -mm  rods 
at  edge,  as  indicated.  Darker  grayscale  values  indicate  areas  where  activity  is 
present  (normalized  to  1  in  this  case);  white  areas  (i.e.,  the  rods)  are  cold. 

face,  at  the  front  face  of  the  collimator,  to  the  center-of-rota- 
tion  (COR).  Images  were  reconstructed  from  projection  data 
using  an  ordered  subsets  expectation  maximization  (OSEM) 
algorithm  [36],  [37]  implemented  with  a  ray-tracing  capability, 
using  8  subsets,  5  iterations,  and  2.5-mm  grid  voxels.  Since 
spatial  resolution  was  not  modeled  within  these  reconstructions, 
there  is  a  limit  to  how  much  spatial  resolution  the  reconstruc¬ 
tions  can  recover,  and  that  limit  will  be  reached  within  a 
relatively  small  number  of  iterations.  The  number  of  subsets 
was  chosen  based  on  previous  evaluations  of  tradeoffs  between 
signal-to-noise  ratio  (SNR)  and  contrast  of  small  lesions  [20], 
[21]. 

A  7.7-cm  diameter  cylinder  containing  a  resolution-fre¬ 
quency  phantom  consisting  of  acrylic  mini-rods  each  2.6  cm 
long  (model  ECT/DLX/MP,  Data  Spectrum  Corp.,  Hills¬ 
borough,  NC)  was  used  as  another  measure  of  the  tomographic 
spatial  resolution  of  the  system  (Fig.  2).  In  each  of  six  sectors, 
the  rods  had  equal  diameters  of  4.7,  3.9,  3.1,  2.3,  1.5,  and 
1.1  mm,  spaced  on  twice  their  diameter.  Water  containing 
15.6  mCi  of  Tc-99m  filled  the  interstitial  spaces  between  the 
rods,  thus,  images  are  of  “cold”  rods.  The  phantom  was  oriented 
such  that  the  rods  were  vertical,  parallel  to  the  axis  of  rotation. 
The  camera  and  phantom  were  both  leveled,  and  the  rods  were 
placed  in  the  central  portion  of  the  field  of  view  (FOV).  The 
ROR  was  fixed  at  4.4  cm.  The  16%  energy  window  yielded 
count  rates  of  ~5.9  kcps.  Simple  circular  tomography  was  used 
for  these  measurements  with  256  projections  over  360°  and  an 
acquisition  time  of  28  s/projection. 


Energy  (keV) 


%  Energy  Resolution  (FWHM/H) 


Fig.  3.  (Top)  Gain  balanced,  summed  energy  spectra  of  all  active  pixels.  As 
per  NEMA-2001  specifications,  the  measured  FWHM  energy  resolution  was 
6.8%.  (Bottom)  Histogram  of  FWHM  Energy  Resolutions  for  all  discrete  pixels. 
Mean  =  6.7%;  Median  =  6.4%;  Min  =  3.5%;  Max  =  24.8%. 


TABLE  I 

Average  Planar  Spatial  Resolution  FWHM  and  FWTM  as  Per  [38] 


Direction  Along  Detector 

Collimator-Source  Horizontal  Vertical 


Distance 

FWHM 

FWTM 

FWHM 

FWTM 

(cm) 

(mm) 

(mm) 

(mm) 

(mm) 

1 

3.4  ±  0.7 

6.3  ±  0.7 

3.7  ±  0.7 

6.6  ±  0.4 

10 

6.7  ±  0.1 

12.2  ±  0.3 

7.3  ±  0.2 

13.0  ±  0.3 

B.  Results 

Of  the  detector’s  5120  pixels,  99.8%  were  routinely  usable. 
The  dead  pixels  were  randomly  dispersed  throughout  the 
imaging  FOV,  with  only  one  bad  4-pixel  cluster  at  the  edge 
of  the  distal  FOV.  The  system  has  a  dead-pixel  correction 
scheme  supplied  by  the  manufacturer,  which  includes  a  linear 
interpolation  of  the  surrounding  pixels’  values.  For  our  imaging 
studies,  the  phantoms  did  not  extend  into  that  portion  of  the 
detector.  One  might  expect  a  similar  dead  pixel  fraction  and 
random  distribution  from  similar  imaging  devices.  Further¬ 
more,  extensive  quality  control  studies  are  performed  before 
each  set  of  experiments;  should  excessive  dead  pixels  appear 
to  be  concentrated  in  any  one  region,  the  modules  could  be 
replaced  within  the  camera.  The  measured  uniformity  of  the 


detector  was  ±4.0%,  given  as  the  overall  standard  deviation  of 
all  pixels  divided  by  the  global  mean. 

The  gain  balanced,  summed  energy  spectra  (Fig.  3)  yielded  an 
intrinsic  energy  resolution  of  6.8%  FWHM  at  140  keV.  Using 
the  same  protocol  for  calculation  of  the  FWHM  (as  per  [38])  for 
each  individual  pixel  resulted  in  a  mean  intrinsic  energy  resolu¬ 
tion  of  6.7%,  with  a  range  of  3.5%  to  24.8%  (Fig.  3,  histogram). 
The  measured  sensitivity  of  the  system  for  the  parallel-beam 
collimator  using  Tc-99m  was  1402  cps/mCi  (37.9  cps/MBq). 

Average  planar  spatial  resolution  FWHM  and  full-width- 
tenth-maximum  (FWTM)  values  as  measured  with  the  capillary 
line  source,  according  to  [38],  for  collimator-source  distances 
of  1  and  10  cm  are  shown  in  Table  I.  Planar  spatial  resolution 
plotted  as  a  function  of  distance  to  the  detector  compared  well 
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Fig.  4.  Planar  spatial  resolution  results  from  capillary  line  source 
measurements  match  calculated  values. 


Fig.  5.  Reconstructed  SPECT  spatial  resolution  results  measured  through  two 
capillary  tubes  located  at  the  ROR  center  and  2-cm  offset.  Data  obtained  from 
10-pixel-wide  profiles  drawn  in  both  the  coronal  and  sagittal  slices  from  the  fifth 
iteration  of  the  reconstructed  data. 

to  the  calculated  system  resolution  (Fig.  4).  Horizontal  and 
vertical  measurements  may  be  different  due  to  the  line  source 
being  ideally  aligned  with  the  collimator  holes  in  one  case, 
and/or  nonalignment  with  the  Cartesian  distribution  of  the 
quantized  CZT  elements.  Furthermore,  the  size  and  shape  of 
the  collimator  holes  are  not  matched  to  the  detector  elements, 
and  the  collimator  holes  are  not  aligned  with  the  Cartesian  dis¬ 
tribution  of  the  square  pixels.  The  system’s  current  unmatched 
collimator  can  potentially  introduce  aliasing  and  sensitivity 
fluctuations  across  an  individual  detector  element.  While  some 
groups  are  investigating  un/matched  collimator  designs  for 
pixellated  detectors  [39]— [41],  it  is  known  that  rotating  the 
unmatched  collimator  with  respect  to  the  detector  pixels  would 
reduce  these  effects  [42]. 

SPECT  spatial  resolution  results  are  shown  in  Fig.  5  for  the 
fifth  iteration  data,  where  10-pixel- wide  profiles  were  drawn 
in  both  the  coronal  and  sagittal  slices,  and  Gaussian  curves 
were  fit  to  the  profiles  to  extract  the  FWHM  information. 


Fig.  6.  Dedicated  emission  mammotomography  system  with  camera  mounted 
on  the  hemispherical  positioning  gantry.  A  filled  anthropomorphic  torso 
phantom  and  pendant,  uncompressed  1060  ml  breast  are  suspended  within  the 
camera’s  FOV. 

Resolution  measurements  are  not  significantly  different  with 
fewer  iterations. 

Reconstructed  results  with  B  subsets,  20  iterations,  and 
2.5-mm  voxels,  from  the  mini-rod  measurements  are  shown  in 
Fig.  2.  The  second  sector  of  rods  (3.9-mm  diameter)  are  clearly 
distinguishable,  as  would  be  expected  with  results  from  Figs.  4, 
and  5  given  the  ROR  of  4.4  cm.  The  3.1 -mm  rods  are  also 
discernible;  though  the  profile  reveals  a  relatively  low  contrast, 
the  rods  can  indeed  be  resolved. 

HI.  Three-Dimensional  Mammotomography 

The  significance  of  any  3-D  orbit  with  a  dedicated  breast 
imaging  system  is  that  it  allows  for  fully  3-D  imaging  of  an  un¬ 
compressed,  pendant  breast  while  avoiding  physical  features  or 
other  hindrances  from  a  patient  during  the  acquisition.  The  to¬ 
mographic  gantry  (Fig.  6)  used  to  position  the  camera’s  central 
ray  at  any  point  in  a  hemisphere  about  the  COR  has  been  pre¬ 
viously  described  [20],  [21].  A  basis  set  of  orbits  was  initially 
characterized  using  a  smaller,  12.8  x  12.8  cm2  FOV  Nal(Tl) 
camera  incorporated  onto  the  dedicated  emission  mammotomo- 
graph  [21],  While  results  showed  significant  improvements  in 
contrasts  and  SNRs  for  breast  lesions  as  compared  to  simple 
vertical-axis-of-rotation  acquisitions  and  uncompressed  planar 
imaging,  the  system  was  limited  by  the  relatively  small  FOV. 
Image  artifacts  resulted  from  a  necessary  camera  shift  to  com¬ 
pensate  for  large,  truncated  breast  volumes. 

The  new  CZT  detector  used  here  has  a  larger  FOV  and  Obvi¬ 
ates  this  camera  shift  of  our  previous  work.  Thus,  in  this  work, 
we  evaluate  a  subset  of  orbits  from  [21]  here  acquired  with  the 
full  FOV  CZT  detector.  In  addition,  we  evaluate  newer  orbits. 
Ideally,  an  acquisition  orbit  should:  1)  meet  Orlov’s  sampling 
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Fig.  7.  (a)  Scale  schematics  of  orbits,  with  1060  ml  breast  modeled.  The 
wire-frame  box  represents  the  16  x  20  cm2  active  detector  FOV,  with  center  on 
the  trajectory  (black  orbit  path)  about  the  breast.  Orbits  are  (Top  to  Bottom): 
Tilted-Parallel-Beam  (TPB);  Circle-Plus-T\vo-Arcs  (CP2A);  Circle-Plus-Arc 
(CPA);  Cloverleaf  (CLOVER);  and  Inverse  Cloverleaf  (invCLOVER)  (see,  also, 
Table  III).  Labeled  in  TPB  schematic  are  the:  center  of  rotation  (COR),  located 
~2  cm  inside  the  breast  along  the  nipple-chest  axis;  direction  of  increased 
polar  camera  tilt  (small  arrow).  Dark  gray  line  from  COR  is  perpendicular  to 
the  detector  face.  Note  that  while  these  implementations  meet  Orlov’s  criteria 
for  sufficient  sampling,  simple  modifications  allow  for  more  viewable  breast 
volume  when  physical  hindrances  preclude  sufficient  sampling.  Next  to  each 
3-D  schematic  are  polar  plots  of  camera  polar  tilt  (<f>)  (plotted  as  a  radius,  with 
90°  at  center  and  0°  at  edge)  versus  azimuthal  angle  (9)  (plotted  around  the 
circle  from  0°  to  360°).  (b)  For  CLOVER  orbit,  example  polar  plot  of  ROR 
contouring  (0  cm  radius  at  center),  which  allows  camera  to  move  close  to  the 
breast,  as  a  function  of  azimuthal  angle  (9). 


criterion  for  sufficient  data  collection  for  accurate  reconstruc¬ 
tion  of  the  activity  distribution  within  the  breast  when  the  ac¬ 
tivity  is  fully  within  the  camera’s  FOV  for  each  projection  angle 
[32];  2)  maximize  the  viewable  breast  volume;  3)  minimize  the 


935  ml  breast  1060  ml  breast 


Fig.  8.  Photographs  of  the  (left)  935  ml  and  (right)  1060  ml  anthropomorphic 
breast  phantoms.  Dimensions  of  both  phantoms  are  given  in  Table  II.  Two  sets 
of  embedded  lesions  were  inserted  into  each  phantom  as  shown. 


TABLE  II 

Dimensions  of  the  Custom  Breast  Phantoms 


Breast 

Dimensions  (in  cm) 

Volume  (ml) 

Superior-Inferior 

Nipple-Chest 

Medial-Lateral 

935 

20 

7.5 

15 

1060 

i _ i* 

12 

18 

background  contamination  from  cardiac  and  hepatic  sources; 
and  4)  contour  the  breast  to  minimize  the  distance  related  reso¬ 
lution  degradations. 

If  used  as  a  secondary  diagnostic  tool,  the  emission  mammo- 
tomography  system  evaluated  here  can  image  an  indeterminate 
lesion  which  may  have  presented  on  an  initial  mammogram.  Or¬ 
bits  can  be  designed  to  contour  the  breast  keeping  the  detector 
position  near  such  a  lesion.  Thus,  for  this  set  of  experiments,  we 
examine  those  orbits  illustrated  in  Fig.  7(a).  Given  the  geometry 
of  a  pendant  breast  and  the  tilts  afforded  by  the  system  to  con¬ 
tour  the  breast,  proximity  of  different  parts  of  the  camera  will 
vary  with  respect  to  the  breast  and  an  azimuthal  axis  of  rotation. 
Orbits  are  designed  to  optimally  contour  the  breast  utilizing  dy¬ 
namic  ROR  control  [example  given  in  Fig.  7(b)],  allowing  at 
most  a  1-cm  separation  between  the  camera’s  center  and  the 
breast  phantom  surface. 

While  this  paper  focuses  on  parallel-beam  imaging  utilizing  a 
focal  spot  at  infinity,  Orlov’s  sampling  criteria  have  been  shown 
to  be  a  special  case  [33]  of  Tuy’s  condition  for  cone-beam  sam¬ 
pling  [43],  which  outlines  criteria  for  cone-beam  type  geome¬ 
tries  that  are  also  possible  with  our  system.  The  flexibility  of 
the  gantry  allows  for  investigation  of  various  sampling  schemes, 
for  a  variety  of  converging  or  even  diverging  collimation  beam 
configurations. 

A.  Methods 

Two  acrylic  lesions,  one  8-mm  and  one  5-mm  inner  diameter 
with  ~l-mm  wall  thickness  ( Data  Spectrum  Corp.)  were  in¬ 
serted  into  a  ~935  ml  custom  fabricated,  fillable  anthropomor¬ 
phic  breast  phantom  ( Radiology  Support  Devices  Inc.,  Newport 
Beach,  CA)  as  shown  in  Fig.  8.  The  8-mm  lesion  was  placed 
centrally  in  the  breast,  while  the  5-mm  lesion  was  placed  proxi- 
mally.  The  custom-made  breast  phantom  is  one  of  several,  each 
with  varying  sizes  and  shapes;  this  particular  model  has  dimen¬ 
sions  given  in  Table  II.  The  dimensions  of  these  phantoms  are 
within  the  ranges  measured  clinically  [44],  Specific  orbit  pa¬ 
rameters  used  for  these  acquisitions,  including  ROR  ranges,  are 
given  in  Table  m.  The  lesion  :breast  radionuclide  concentra¬ 
tion  ratio  was  9.8:1,  with  an  absolute  lesion  activity  concen¬ 
tration  of  12.1  /iCi/ml.  The  initial  scan  was  acquired  for  20 
min,  with  subsequent  acquisitions  adjusted  to  compensate  for 
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TABLE  III 

Orbit  Parameters  Used  for  Each  Breast  Phantom  Acquisition, 
With  All  Orbits  Acquired  Over  a  360°  Azimuthal  Range  (0) 


Orbit  Name 

Acronym 

< p  Range 
(min-max) 
(degrees) 

ROR  Range 
(min-max) 
(cm) 

Number 

of 

Ptjs 

935  ml  Breast,  COR 

ocated  at  nipple 

Tilted-Parallel-Beam 

TPB 

45° 

1.9  -  3.7 

128 

Circle-Plus-Two- Arcs 

CP2A 

0-45° 

1.9  -  9.3 

162 

Circle-Plus- Arc 

CPA 

0  -  45° 

2.4  -  9.3 

145 

Cloverleaf 

CLOVER 

0-45° 

2.0  -  9.3 

268 

Inverse  Cloverleaf 

invCLOVER 

0  -  45° 

2.0  -  9.3 

268 

1060  ml  Breast  +  Filled  Torso,  COR 

1  located  ~2  cm  inside  nipple 

Tilted-Parallel-Beam 

TPB 

45° 

3.2  -  4.4 

128 

Circle-Plus-Arc 

CPA 

15  -  45° 

3.8  -  7.1 

139 

Cloverleaf 

CLOVER 

15  -  45° 

3.6  -  7.2 

256 

Inverse  Cloverleaf 

invCLOVER 

15  -  45° 

3.6  -7.1 

256 

TABLE  IV 

Activity  Concentration  Ratios  for  Each  of  the  Filled  Organs 
in  the  Anthropomorphic  Torso  and  Breast  Phantoms 


Organ 

Volume 

(ml) 

Concentration 

(/rCi/ml) 

Ratio 

Breast 

1060 

0.33 

1.0 

Torso 

7810 

0.33 

1,0 

Liver 

1010 

4.2 

12.7 

Heart  (outer  wall) 

290 

4.0 

12.1 

9.5  mm  diam.  lesion 

0.43 

2.0 

6.1 

8.0  mm  diam.  lesion 

0.26 

2.0 

6.1 

radioactive  decay.  As  in  Section  II,  a  ±8%  symmetric  energy 
window  about  the  140-keV  photopeak  was  used,  yielding  count 
rates  of  ~0.9  kcps.  The  breast  phantom  was  attached  to  the  an¬ 
thropomorphic  torso  phantom  to  simulate  complications  in  pa¬ 
tient  positioning,  but  the  torso  was  not  filled  with  radioactivity. 

To  more  nearly  simulate  a  clinical  study  and,  thus,  match  the 
biodistribution  of  activity  in  a  patient,  an  anthropomorphic  torso 
phantom  containing  a  heart,  liver,  and  lungs  ( Radiology  Support 
Devices  Inc)  was  filled  according  to  compiled,  published  data 
[3],  [4],  [45]— [47].  The  ~1060  ml  breast  phantom,  torso  and  or¬ 
gans  were  filled  according  to  the  concentration  ratios  given  in 
Table  IV,  assuming  an  initial  total  dose  of  25  mCi  (925  MBq), 
with  the  heart  uptake  of  1.2%  after  5  min  [47],  Two  lesions, 
one  9.5-mm  and  one  8.0-mm  inner  diameter  with  ~l-mm  wall 
thickness  (Data  Spectrum  Corp .)  were  also  filled  according  to 
the  concentration  ratios  given  in  Table  IV.  The  9.5-mm  lesion 
was  placed  centrally  in  the  breast,  while  the  8-mm  lesion  was 
placed  in  the  upper  outer  quadrant  of  the  breast  (Fig.  8).  Total 
image  time  was  adjusted  to  compensate  for  radioactive  decay, 
with  the  initial  scan  acquired  for  20  min.  Because  of  the  ex¬ 
pected  increase  in  secondary  (scattered)  contamination  resulting 
from  the  filled  torso,  a  tighter  energy  window  was  used  to  reject 
scatter.  Thus,  in  these  breast+filled  torso  acquisitions,  a  sym¬ 
metric  ±4%  window  about  the  140-keV  photopeak  of  Tc-99m 
was  used,  providing  count  rates  of  ~0.2  kcps. 

For  the  3-D  camera  trajectories  implemented  in  these 
breast+filled  torso  experiments,  the  filled  torso  distal  to  the 
breast  introduces  background  activity  which  is  not  uniformly 
in  the  FOV  for  each  projection  angle.  Hence,  Orlov’s  criteria 
is  not  strictly  met.  Thus,  the  primary  design  goal  for  the 
acquisition  orbits  shifts  from  meeting  Orlov’s  criterion  to  max¬ 
imizing  the  viewable  breast  volume.  In  order  to  increase  the 
viewable  volume,  the  COR  was  located  as  far  inside  the  breast 


as  practicable  to  image  near  the  chest  wall.  As  a  result,  these 
orbits  have  a  minimum  polar  tilt  of  15°,  with  the  exception  of 
tilted-parallel-beam  (TPB)  which  maintains  a  larger,  45°  polar 
tilt  (Table  III). 

Images  were  again  reconstructed  from  projection  data  using 
the  OSEM  algorithm  with  ray-tracing  capability,  with  8  subsets, 

5  iterations,  calculated  attenuation  correction  assuming  a  uni¬ 
form  emission  volume,  and  2.5-mm  grid  voxels.  A  3-D  Hann 
filter  with  a  cutoff  of  0.7  times  the  Nyquist  frequency  was  ap¬ 
plied  to  the  reconstructed  data. 

Regions  of  interest  (ROIs)  were  drawn  both  inside  the  le¬ 
sions  and  circumferentially  about  each  lesion  in  the  breast  back¬ 
ground.  For  lesion  sizes  of  9.5-,  8-,  and  5-mm-diameter,  ROIs 
contained  12,  9,  and  4  pixels,  respectively.  In  all  cases,  24  sep¬ 
arate  background  regions  were  drawn  about  each  lesion,  8  each 
in  the  slice  containing  the  lesion’s  centroid,  in  the  2.5-mm  slice 
above,  and  in  the  2.5-mm  slice  below.  The  size  of  the  back¬ 
ground  ROIs  drawn  about  each  lesion  corresponded  to  the  size 
of  the  ROI  drawn  within  that  lesion  (e.g.,  for  the  9.5-mm  lesion, 
24  background  regions  totaled  288  pixels,  whereas  for  the  8-mm 
lesion,  24  background  regions  totaled  216  pixels). 

The  SNR  was  measured  as  the  difference  in  mean 
pixel  values  of  the  lesion  and  breast  background,  divided 
by  the  standard  deviation  of  the  uniform  background 
((ROIicsion  -  ROIbreast )  / ^breast )  •  Lesion  contrast  is  de¬ 
fined  as  the  signal  difference  divided  by  the  background  signal 

(  (  R  O  lies-'nn  ROIbreast  )/IIDIbrcast)- 

B.  Results 

Fig.  9  shows  reconstructed  results  from  the  935  ml  breast- 
only  measurements,  where  both  lesions  are  located  in  the  same 
transaxial  plane.  Without  additional  torso  contamination,  and 
with  a  uniform  breast  background,  both  the  centrally  located 
8  mm  and  the  proximally  located  5-mm-diameter  lesions  are 
clearly  visible.  Table  V  gives  results  from  calculations  of  SNR 
and  contrast  for  both  lesions  at  the  second  OSEM  iteration, 
smoothed  with  a  Hann  filter  (which  was  the  optimal  tradeoff 
between  gains  in  contrast  and  losses  in  SNR  (see  [20],  [21])). 

Fig.  10(a)  shows  reconstructed  results  from  the  1060  ml 
breast  and  filled  anthropomorphic  torso  phantom  measure¬ 
ments,  for  the  centrally  located  9.5-mm-diameter  lesion. 
Fig.  10(b)  shows  an  example  of  reconstructed  results  for  the 
8-mm  lesion.  Table  V  lists  results  from  calculations  of  SNR 
and  contrast  for  both  the  9.5-mm  and  8-mm  lesions,  again  at 
the  second  OSEM  iteration,  smoothed  with  a  Hann  filter. 

C.  Discussion 

The  two  breast  phantoms  were  chosen  for  these  mammoto- 
mography  studies  because  of  their  similar  size  (volume),  but 
vastly  different  shape.  The  breast’s  shape  was  purposely  and 
deliberately  altered  to  examine  its  effect  on  the  newer  orbits 
evaluated  within  this  work.  With  a  different  breast  shape,  var¬ 
ious  complications  in  positioning  arise  (including  ROR  control 
to  contour  the  breast)  without  significantly  changing  the  amount 
of  scatter,  total  activity,  or  attenuation  that  one  might  expect  with 
vastly  different  sized  breasts. 

In  studies  with  the  935  ml  breast,  high-count,  low-noise  im¬ 
ages  result  (Fig.  9).  This  is  not  surprising  given  the  fairly  large 
energy  window  (±8%)  and  large  absolute  activity.  However, 
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Fig.  9.  OSEM  reconstructed  data  (8  subsets,  second  iteration  smoothed  with 
Hann  filter  shown)  from  the  935  ml  breast  data,  for  each  labeled  orbit,  with  a 
9.8:1  lesiontbreast  radionuclide  concentration  ratio.  Solid  line  in  transaxial  slice 
and  profile  is  through  the  8-mm  lesion,  and  the  dashed  line  is  through  the  5-mm 
lesion.  Note  that  incomplete  sampling  artifacts  from  reconstructed  data  seen 
with  the  TPB  orbit  sagittal  views  are  reduced  with  the  more  complex  orbits  that 
improve  sampling,  similar  to  results  shown  in  [20]  and  [21], 

this  ideal  case  does  demonstrate  the  ability  of  the  system  to 
clearly  resolve  a  5-mm  lesion.  Furthermore,  both  subcentimeter 
diameter  lesions  were  resolved  with  each  improved  3-D  sam¬ 
pling  scheme,  afforded  by  the  flexibility  of  the  hemispherical 
positioning  gantry,  suggesting  the  viability  of  these  trajectories 
for  clinical  studies.  With  the  use  of  a  radio-opaque  bed  to  shield 
the  system  from  additional  body  background,  we  are  currently 
investigating  whether  the  image  quality,  given  enough  statistics, 
could  ideally  resemble  that  shown  in  these  isolated  breast  mea¬ 
surements  [48], 

The  TPB  orbit  (which  does  not  strictly  meet  Orlov’s  suffi¬ 
ciency  condition)  yields  the  lowest  SNR  value  for  the  5-mm 
lesion.  Further,  the  reconstructed  images  are  riddled  with  arti¬ 
facts,  including  an  elongated  nipple-chest  axis,  similar  to  our 


TABLE  V 

Measured  SNRs  and  Contrasts,  at  2nd  OSEM  Iteration, 
Smoothed  With  a  Hann  Filter 


9.5  mm 

Lesion  Size 

8  mm 

5  mm 

(0.43  ml) 

(0.26  ml) 

(0.06  ml) 

Orbit 

SNR  Contrast 

SNR  Contrast 

SNR  Contrast 

Ideal  low-noise  case  with  high  uptake: 

935  ml  Breast  Only 

10:1  activity  concentration  ratio,  12.1  pCilml  in  lesions,  16%  EWin 


central 

distal 

TPB 

- 

- 

19.3 

2.0 

2.7 

0.2 

CP2A 

- 

- 

26.6 

1.9 

4.9 

0.2 

CPA 

- 

— 

17.3 

1.4 

3.9 

0.2 

CLOVER 

_ 

- 

23.2 

1.4 

7.9 

0.3 

invCLOVER 

- 

17.1 

1.5 

6.9 

0.3 

More  realistic  case  simulating  patient  biology: 

1060  ml  Breast  +  Filled  Torso, 

6.1:1:12  activity  concentration  ratio,  2.0  ptCi/ml  in  lesions,  8%  EWin 


central 

distal 

TPB 

4.6 

0.8 

3.4 

0.5 

- 

- 

CPA 

4.8 

1.0 

5.5 

0.6 

- 

- 

CLOVER 

3.7 

0.8 

2.0 

0.3 

- 

- 

invCLOVER 

6.3 

1.2 

1.2 

0.2 

- 

- 

previous  results  [20],  [21].  With  the  more  complex  orbits  that 
improve  volumetric  sampling,  the  breast  shape  is  almost  fully 
recovered.  Results  also  indicate  that  the  orbit  most  optimal  for 
specific  lesion  locations  may  change;  whether  this  result  is  more 
ascribable  to  the  placement  of  the  lesion  or  the  size  of  the  lesion 
itself  is  being  further  investigated. 

The  viewable  breast  volume  is  limited  for  the  935  ml  breast 
case  along  the  nipple-chest  axis.  Because  the  design  criteria  in¬ 
cluded  meeting  Orlov’s  sampling  criteria,  the  COR  was  located 
at  the  nipple  for  this  particular  breast  shape  and  size.  Thus,  with 
the  lesions  optimally  placed  near  the  nipple  for  this  set  of  mea¬ 
surements,  resulting  images  provided  for  high  contrast  resolu¬ 
tion.  However,  for  lesions  located  near  the  chest  wall,  relatively 
far  from  the  detector,  a  COR  located  at  the  nipple  may  pre¬ 
clude  their  visualization  even  given  improved  sampling.  Conse¬ 
quently,  design  criteria  for  any  orbit  require  a  tradeoff  between 
improved  sampling  and  viewable  breast  volume. 

This  subsequent  optimization  of  orbit  design  parameters  was 
examined  in  the  1060  ml  breast  +  filled  torso  experiments.  The 
design  goal  for  the  orbit  implementations  for  those  measure¬ 
ments  was  to  move  the  COR  as  far  into  the  breast  as  practi¬ 
cable.  Accordingly,  the  minimum  tilt  was  increased  from  0°  to 
15°  so  that  the  COR  could  be  moved  into  the  breast  ~2  cm. 
This,  in  turn,  increased  the  viewable  breast  volume.  Indeed,  the 
entire  12-cm  length  of  the  nipple-chest  wall  axis  can  be  seen, 
as  well  as  some  additional  depth  into  the  chest  wall  with  only 
a  small  missing  cone  of  data  behind  the  chest  wall  along  the 
nipple-chest  axis. 

With  more  realistic  uptake  ratios  in  the  body  and  absolute 
activities  of  the  organs,  the  9.5-mm  centrally  located  lesion 
(~6  cm  from  the  nipple)  was  discernible  in  all  images.  Again, 
results  from  the  TPB  orbit  in  Fig.  10(a)  demonstrate  the  familiar 
elongation  of  the  nipple-chest  axis.  Additionally,  increased  ac¬ 
tivity  from  the  background,  specifically  the  liver,  is  also  visible 
[Fig.  10(a),  white  arrow].  Interestingly,  the  remaining  orbits, 
though  also  incompletely  sampling  the  breast  according  to 
Orlov’s  condition,  appear  to  more  nearly  recover  the  breast’s 
shape.  This  suggests  that  while  the  mathematical  condition 
for  accurate  reconstruction  of  the  activity  distribution  within 
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Fig.  10.  OSEM  reconstructed  data  (8  subsets,  second  iteration  smoothed 
with  Hann  filter  shown)  from  the  1060  ml  breast  +  filled  torso  data,  for 
each  labeled  orbit  with  an  incomplete  sampling  scheme,  with  a  6.1:1:12 
lesion:breast-and-torso:liver-and-heart  radionuclide  concentration  ratio, 
(a)  Solid  line  is  drawn  through  the  9.5-mm  lesion.  Note  the  increased 
background  contamination  with  TPB  (white  arrow)  as  compared  to  the  other 
orbits  which  do  not  maintain  the  large  detector  polar  tilts,  thus  limiting  more 
direct  views  of  the  heart  and  liver.  Even  with  these  incompletely  sampled 
orbits,  and  more  realistic  absolute  radioactivities  within  the  phantom,  the 
centrally  located  9.5-mm  lesion  is  clearly  visible.  Dashed  arrow  in  CPA  points 
to  cold  lesion  holder,  which  can  also  be  seen  with  the  other  orbits  (see,  also, 
InvCLOVER)  and  in  Fig.  8.  (b)  Views  of  the  8-mm  lesion  at  different  transaxial 
and  sagittal  locations  than  the  larger  lesion. 


the  breast  was  not  met,  in  practice,  the  sampling  requirement 
may  be  less  strict  than  initially  anticipated;  this  is  an  area  that 
warrants  further  investigation. 

Incomplete  sampling  may,  however,  produce  false  positives 
in  the  image  [Fig.  10(a),  CLOVER];  this  may  also  be  due  to 


other  reconstruction  artifacts  since  these  misleading  signals  are 
actually  manifested  as  streaks  and  appear  in  multiple  contiguous 
planes.  The  more  peripheral  areas  of  the  breast  along  the  chest 
wall  may  suffer  from  less  sufficient  sampling;  the  uniformity 
of  the  sampling  scheme  itself  is  currently  being  investigated  in 
addition  to  completeness  issues. 

For  this  work,  there  was  no  patient  (torso)  shielding  whatso¬ 
ever.  Given  a  radio-opaque  support  bed,  however,  background 
activity  should  be  significantly  reduced,  as  the  heart  and  liver 
should  be  removed  from  a  significant  number  of  detector  van¬ 
tages.  The  shielded  bed  may  also  affect  the  allowed  camera 
trajectories.  Our  group  is  in  the  process  of  implementing  and 
investigating  such  a  patient  bed  and  its  effects  (and  expected 
improvements)  on  image  quality  [48]. 


IV.  Conclusion 

Dedicated  emission  mammotomography  with  a  CZT  imaging 
detector  that  could  image  large  breasts  was  implemented.  Ini¬ 
tial  characterization  and  preclinical  evaluation  studies  were 
performed.  The  intrinsic  characteristics  of  the  medium-FOV, 
semiconductor  imaging  detector  were  evaluated,  as  was  the 
overall  suitability  of  the  device  for  dedicated  SPECT  imaging 
of  the  breast.  The  detector  was  found  to  have  a  6.8%  FWHM 
at  140  keV,  and  was  able  to  clearly  discern  3.1-mm-diameter 
mini-cold  rods  in  reconstructed  images.  A  key  benefit  provided 
by  a  dedicated,  compact  SPECT  breast  imaging  system  is 
its  ability  to  acquire  fully  3-D  molecular  data  about  an  un¬ 
compressed,  pendant  breast  and  possibly  anterior  chest  wall. 
As  a  potential  secondary  diagnostic  tool,  the  3-D  sampling 
scheme  can  be  tailored  to  an  individual  patient,  especially  if 
mammographic  studies  can  provide  a  priori  information  about 
a  suspected  lesion.  Though  it  appears  no  single  orbit  optimizes 
image  quality  for  lesions  at  all  placements  within  the  breast, 
the  phantom  studies  shown  here  demonstrate  the  capability  of 
the  system  to  easily  adapt  the  acquisition  scheme  for  imaging 
of  various  breast  shapes  and  sizes.  Additionally,  there  are 
no  outstanding  trajectories  that  consistently  yield  optimized 
quantitative  lesion  imaging  parameters.  Qualitatively  however, 
while  imaging  isolated  breasts  which  are  completely  in  the 
field-of-view  for  all  vantages  yields  consistent  reconstructed 
images,  imaging  breasts  with  realistic  torso  backgrounds 
(out-of-field  activity)  substantially  alters  image  characteristics 
and  breast  morphology.  Orbits  that  improve  sampling  yield 
more  realistic  morphology  as  well  as  more  accurate  overall 
activity  distributions  throughout  the  imaged  volumes.  Fur¬ 
thermore,  results  from  breast  phantom  measurements  reveal 
the  ability  of  this  approach  to  visualize  lesions  <1  cm  in 
diameter,  and  studies  are  ongoing  to  evaluate  this  ability  with 
vastly  different  sized  breast  phantoms  [48].  The  strength  of  this 
molecular  imaging  approach  for  breast  imaging  can  be  further 
enhanced  with  additional  structural  information  provided  by 
transmission  cone-beam  mammotomography.  Initial  devel¬ 
opment  of  a  prototype  X-ray  computed  mammotomograph 
for  such  a  purpose  is  underway  [30],  [31].  The  3-D  emission 
mammotomography  system  evaluated  within  this  work,  incor¬ 
porating  a  compact  CZT  detector,  shows  a  promising,  clinically 
relevant  technique  for  molecular  imaging  of  the  whole  breast. 
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